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Abstract 

The observational progress in the y-ray astronomy in the last few years has led to 
the discovery of more than a thousand sources at GeV energies and more than a 
hundred sources at TeV energies. A few different classes of compact objects in 
the Galaxy have been established. They show many unexpected features at high 
energies the physics of which remains mainly unknown. At present it is clear that 
detailed investigation of these new phenomena can be performed only with the 
technical equipment which offer an order of magnitude better sensitivity, and a 
few times better energy, angular and time resolution in the broad energy range 
staring from a few tens of GeV up to a few hundreds TeV. Such facilities can 
be realized by the next generation of instruments such as the planned Cherenkov 
Telescope Array (CTA). 

The aim of this report is to summarize up to date observational results on the 
compact galactic sources in the GeV-TeV y-ray energy range, discuss their theo- 
retical implications, and indicate which hypothesis considered at present might be 
verified with the next generation of telescopes. We point out which of the obser- 
vational features of the y-ray sources are important to investigate with special care 
with the planned CTA in order to put a new light on physical processes involved. 
Their knowledge should finally allow us to answer the question on the origin of 
energetic particles in our Galaxy. 
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1. Introduction 



It is widely believed that relativistic, charged particles are accelerated in the 
Galaxy up to at least ~10 17 eV. Such particles are confined within the Galactic 
disk since their Larmor radius is smaller than the thickness of the Galactic disk. 
However, which types of the observed sources are responsible for the acceleration 
processes (maybe a few ?) is at present not quite clear. The best candidates seem to 
be sources containing compact objects such as rotating neutron stars (NSs), solar 
mass black holes, massive stars, or maybe also White Dwarfs (WDs). They can 
either form separate systems which energize their surrounding (e.g. Pulsar Wind 
Nebulae) or occur close to concentrations of matter within the Galaxy (Globular 
and Open Clusters) or be companions of normal stars (binary systems). 

During the last few years our knowledge on the compact sources of cosmic 
accelerators in the Galaxy has experienced a major breakthrough. Thanks to 
the development of a new satellite and on-ground telescopes, a few hundred of 
discrete y-ray sources at energies above 100 MeV have been discovered in the 
Galactic plane by the Large Area Telescope on board of Fermi satellite (Fermi- 
LAT) Q, and a few tens of sources above ~100 GeV have been discovered 
by the Cherenkov telescopes (see for the reviews, e.g. [|182[ 1471 111911 . Some of 
these sources have showed unexpected emission features which put new light on 
the particle acceleration and radiation processes occurring in them. 

One of the purposes of this paper is to present up-to-date knowledge (both 
observational and theoretical) on the emission features of at present very well es- 
tablished y-ray sources such as: different types of pulsars, Pulsar Wind Nebulae 
(PWNe), Globular Clusters (GCs) and binary systems containing either two mas- 
sive stars, a WD, a pulsar, or a solar mass black hole. We also mention other 
potential compact y-ray sources in the Galaxy such as: isolated massive stars and 
WDs, or accreting NSs and WDs. They are expected to emit y-rays based on the 
theoretical analysis of the conditions within them which seem to be suitable for 
acceleration of particles and subsequent production of high energy radiation. We 
will consider how our knowledge of high energy processes in the Galactic sources 
can profit from the construction of the next generation telescopes able to observe 
y-rays with energies above ~30 GeV (such as the Cherenkov Telescope Array 
- CTA, the High Altitude Water Cherenkov Experiment - HAWC, or the Water 
Cherenkov Detector Array (WCDA) of the Large Altitude Air Shower Observa- 
tory - LHAASO). For example, CTA is planning to reach 50 hr sensitivity (5 a 
detection) on the level of 1% Crab flux units (CU) at ~30 GeV and above ~30 
TeV, and 0. 1% CU in the range from 100 GeV to 3 TeV H- This will be about an 
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order of magnitude better than achieved with the present telescope systems. The 
angular resolution of CTA is expected to be at least a factor of 2-3 better than that 
of the presently operating HESS Array. On the other hand, HAWC will provide 
all sky monitoring with 15 times better sensitivity than the Milagro Observatory 
(see HAWC Gamma-ray Observatory web page: http://hawc.umd.edu). It is ex- 
pected that instruments with such sensitivities will finally allow us: to answer the 
question about the acceleration sites of particles with energies at most ~10 15 eV, 
to conclude about their acceleration mechanism, and provide constraints on the 
physics of these high energy sources. These new instruments are also expected to 
discover new types of sources, maybe some of them mentioned in this review. 



2. Isolated Rotating Neutron Stars 

Isolated rotating neutron stars (pulsars) produce winds and radiation that pre- 
vents accretion of matter onto their surfaces. They are at present one of the best 
known and defined astrophysical objects observed in the broad energy range, from 
radio up to ~100 GeV. Their parameters can be precisely measured (period, pe- 
riod derivative) making it possible to estimate their basic physical quantities in 
terms of the rotating dipole model such as: the age, the surface magnetic field 
strength, the rotational energy loss rate. These characteristics allow us at present 
to distinguish three basic types of pulsars: (1) classical pulsars (CPs) with the sur- 
face magnetic field close to a few ~10 12 G and rotational periods between several 
milliseconds and a few seconds; (2) millisecond pulsars (MSPs) with the surface 
magnetic fields close to a few 10 8 G and periods of a few milliseconds; and (3) 
anomalous X-ray pulsars (AXPs) and soft y-ray repeaters (SGRs) interpreted in 
terms of the magnetar model in which a slowly rotating neutron star (period of a 
few to several seconds) has a superstrong surface magnetic field of the order of 
~10 14 G. These specific classes of pulsars are not well separated. Objects with 
intermediate parameters are also observed. Their existence can tell us a lot about 
the origin of specific classes, especially in the case of the nature of AXPs and 
SGRs since their magnetar interpretation has sometimes been questioned. Recent 
discoveries of the pulsed y-ray emission at energies >100 MeV from the high 
field magnetic radio p ulsar s, such as PSR B 1509-28 [1 1], the pulsar in CTA 1 



and PSR Jl 1 19-6127 |238|], with the features similar to CPs have strengthened the 
magnetar hypothesis. 

Due to different parameters of these types of pulsars we consider below their 
possible high energy processes separately. For a recent review of the basic prop- 



erties of these pulsars see [|209ll 
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2.1. Classical pulsars 

With the launch of the Fermi and AGILE satellites, the number of known 
classical y-ray pulsars has increased by more than an order of magnitude from 
6 (discovered by EGRET telescope) to almost ~100 [|12|l within only 2-3 years. 
The excellent sensitivity and timing of the Fermi-LAT allowed this collaboration 
to detect several pulsars in blind periodicity searches [4] for the first time. 

2.1.1. New results on old y-ray pulsars 

Studies of the strongest pulsars (detected already by EGRET) showed a va- 
riety of interesting details. For example the two peak structure observed in the 



Crab pulsar in the radio and y-ray s are slightly shifted by ~0.3 ms ||13(|. The Crab 
pulsar light curve shows interpulse emission which has clearly harder spectrum 
than that coming from the peak regions. In the case of the Vela pulsar, the y-ray 
baseline emission extends through 80% of the li ght curve, showing an intriguing 
third peak, more evident at higher energies [R jilll. The first peak in the light 
curve of Vela dominates at low energies, but it disappears completely at the high- 
est energies (above a few GeV). A fit to the whole pulsed y-ray spectrum by an 
exponential function does not reproduce correctly the data at the highest detected 
energies. The third famous pulsar, Geminga, shows many similarities to the pre- 
vious two pulsars, although up to now no radio emission has been observed from 
it. Geminga shows y-ray emission at all phases which is consistent with y-ray 
production at regions close to the light cylinder [15]. Another interesting pulsar 
with quite strong surface magnetic field, PSR B 1509-58, has been originally de- 
tected only by the COMPTEL telescope 1120011 . Now, the pulsed emission is also 



observed by the AGILE and Fermi-LAT at energies below ~1 GeV II235L I11I1. The 
y-ray spectrum shows a sharp break at a few tens of MeV, a feature not observed 
up to now in any other pulsar. 

The whole population of pulsars detected by the Fermi-LAT shows a variety 
of light curves. These light curves can be divided into three main groups: those 
with two peaks well separated (by ~0.4-0.5), with two close peaks (separated by 



only ~0.2) and those with a single broad peak [12]. The differential spectra are 
well described by the power law functions with spectral indices close to 1 .5 and 
exponential cut-offs at a few GeV. At present it is difficult to conclude what is the 
conversion efficiency of the rotational energy loss rate to the y-ray power due to 
the unknown geometry of the y-ray emission and uncertain estimates of distances 
to the pulsars. 
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Crab Pulsar, P1 + P2 
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Figure 1 : The pulsed y-ray spectrum of the Crab pulsar (both peaks together) observed with the 
Fermi, MAGIC and VERITAS telescopes. Specific measurements are described in the legend 
within the figure (from H58I1 . reproduced by permission of the A&A). 



2.1.2, Unexpected detections of the Crab above several GeV 

Pulsars have become interesting targets for the presently operating Cherenkov 
telescopes even before launching the Fermi-LAT detector. For the first time, the 
Crab pulsar was detected at energies above ~25 GeV by the MAGIC Collabora- 
tion, using the special pulsar trigger [16011 . Some evidence of the signal at energies 
above ~60 GeV has been also reported from this pulsar [jHH. The lack of a sharp 
cut-off in the pulsed Crab spectrum at a few GeV has been also confirmed by the 
Fermi-LAT observations up to a few tens of GeV uM- Further observations with 
the MAGIC and VERITAS telescopes allowed the grou ps t o m easur e the pulsed 
y-ray spectrum from the Crab pulsar up to ~400 GeV [62, 57, l260ll . This total 
pulsed y-ray spectrum is well described by a simple power law extending from the 
last Fermi-LAT point (see Fig. 1). The y-ray pulses become narrower for larger 
energies. The ratio of the height of the first pulse (PI) to that of the second pulse 
(P2) in the Crab decreases with energy. 

No other pulsed emission above a few tens of GeV has been reported up to by 
the modern telescope groups (e.g. [|44l l52l |67J. The pulsar spectra measured by 
Fermi-LAT seem to be consistent with exponential cut-offs at a few GeV [12]. 
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2.1.3. Recent observations versus pulsar models 

The y-ray production in pulsars is interpreted in terms of a few general models 
locating the acceleration and y-ray e mission reg ions in different parts of the in- 
ner magneto sphere. Polar cap models 141] assume that acceleration occurs 
close to the polar regions of the NS surface, y-rays are then produced by curva- 
ture radiation ||l4l|| or resonant/non-reso nant Inverse Compton Scattering (ICS) 
of thermal radiation from the NS surface 1129 ill . The spectra are expected to show 
super-exponential cut-off at a few GeV due to the y-ray absorption in a strong 
magnetic field. In a more recent version of this model, the so called pair starved 
polar cap (e.g. Q297[|), particles are accelerated to large distances from the NS sur- 
face since the electric field is not saturated by inefficient production of e ± pairs in 
the cascades. In such a case, y-rays can be produced at larger distances from the 
NS surface, not sufferring strong absorption in the magnetic field. Polar cap model 
can also be extended to the region close to the last open magnetic field lines up to 
the light cylinder where the curvature of magnetic field lines prevents production 
of the next generation of e ± pairs along them flVOL 1 1 VSh . A technical realization of 
this model has been successfully applied to a detailed modelling of a num ber o f 
pulsar y-ray light curves and spectra (the so-called two-pole caustic model 11570. 
In this model the y-ray light curves are explained as the contributions from the 
regions of the two magnetic poles. 

The acceleration of particles and y-ray production can also occur in the outer 
pulsar magnetosphere in the region limited by th e null su rface and the light cylin- 
der, according to the so called outer gap model fll 30L Il3ll I251N . The two versions 
of the model differ in the origin of the e ± pairs which limit the extension of the 
outer gap. In the Cheng et al. model, e pairs are produced by the primary TeV 
y-rays which collide with the synchrotron radiation from the outer gap. In the Ro- 
mani model, the e ± pairs are produced in collisions of primary curvature y-rays 
with thermal X-rays from the NS surface. In fact, the outer gap models predicted 
the existence of the low level, pulsed hard y-ray spectrum extending up to TeV en- 
ergies as a result of the comptonization of the synchrotron radiation by primary e^ 
pairs. At present it is not clear whether this component might be responsible for 
the sub-TeV pulsed emission from the Crab pulsar. Detailed calculations of the ex- 
pected y-ray spectra and the light curves in this model ar e complicated since they 
require a three-dimensional treatment (see e.g. 11278 , Depending on the 

parameters of the pulsars, y-ray spectra are formed in the outer gap model either 
by synchrotron an d curvature processes or, in the case of very young pulsars, also 
by the IC process IU32L l277ll . The above considered models are based on simple 
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pulsar magneto sphere structure (approximated by the vacuum dipole magnetic 
model). Recent numerical studies of rotating force-free pulsar magnetospheres 
enable one determine a magnetic field structure in the outer magnetosphere and, 
as a result, a more detailed calculations of y-ray emission patterns (the so-called 



separatrix Layer model [|75|]). Based on these calculations, it is concluded that the 
two-pole caustic and the outer gap models are not able to reproduce the observed 
pulse structures in the y-ray pulsar light curves. 

Recent y-ray observations made for the first time possible to put strong con- 
straints on the y-ray emission region within the pulsar magnetosphere, locating it 
at least a few stellar radii from the neutron star surface (closer in, the gamma-rays 
would be absorbed in the magnetic field and converted to e ± pairs [78]). This actu- 
ally eliminated the popular scenario for the y-ray origin close to the stellar surface 
by the curvature radiation mechanism, the so-called polar cap model. However, 
these observations are still consistent with the pulsar models in which most of the 
y-rays are produced far away from the neutron star surface, i.e. the slot gap and 
the outer gap models. 

In the context of the recent discovery of the sub-TeV emission from the Crab 
pulsar, a production of y-rays by the secondary cascade e pairs in the IC scat- 
tering of the s ynchrotron UV-X-ray radiation from the gap was considered by 



Lyutikov 1121211 and the detailed calculations along such model were performed by 
Hirotani IU84I1 . see also the comparison with the observations in Q5TD - 

The pulsed sub-TeV component might also appear as a result of comptoniza- 
tion of the magneto spheric soft radiation by e ± pairs in the relativistic pulsar wind 
above the light cylinder [108]. By tuning the distance of the emission region from 
the light cylinder and the Lorentz factor of the wind the emission was obtained in 
the sub-TeV region. The present observations allow us to put constraints on the 
wind formation region in the Crab nebula based on the last model. The pulsed 
y-ray emission was also propos ed to originate in the pulsar wind region in terms 
of the striped wind model II 196U24 ill . The peaks of y-ray emission appear as a re- 
sult of comptonization of the background radiation by electrons moving radially 
outside the light cylinder with the Doppler factor of the striped wind. 

2.1.4. Interesting problems to be studied with CTA 

The features of the Crab pulsar spectrum at sub-TeV energies were not pre- 
dicted by any detailed model and appeared as a surprise for the pulsar community. 
Present observations suggest that the pulsed y-ray spectrum from the Crab pul- 
sar show a smooth transition between the GeV to sub-TeV energy range (above a 
break at a few GeV). The question appears whether the sub-TeV pulsed emission 
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means a new component in the spectrum (another emission process as expected 
in some models?) or it is a continuation of the spectrum from the GeV energy 
range? Detailed observations of the transition between these two components 
with the next generation telescope systems (e.g. CTA) should solve this prob- 
lem (i.e. whether two distinct components are present in the sub-TeV spectrum 
or where is the transition between these components?). It is also very important 
to establish the end of the TeV component in the pulsed Crab spectrum in order 
to constrain the maximum energies of emitting e ± pairs and conclude about the 
exact mechanism of the y-ray production. Are these y-ray produced in the scat- 
tering process of the soft synchrotron radiation from the outer gap or the thermal 
component from the NS surface or maybe in the curvature process of particles 
accelerated close to the light cylinder? 

Other pulsars discovered by the Fermi-LAT seem to have y-ray spectra consis- 
tent with the exponential cut-offs as expected in the curvature radiation process. 
However, they were not observed with the Cherenkov telescopes up to now. One 
of the main tasks of CTA should be an extensive investigation of the large popu- 
lation of the y-ray pulsars in order to determine which of them (maybe all?) show 
some level of sub-TeV emission? It is supposed, in terms of the outer gap model, 
that only very young pulsars with the Crab-like magnetospheres (and two well 
separated peaks) can produce sub-TeV emission. Observation of other classes of 
pulsars (Vela-type, Geminga-type or single pulsed) with CTA should significantly 
constrain the emission models. 

CTA should also be able to investigate in detail additional emission features 
(third peaks, shifts in localizations of peaks observed in different energies) in the 
pulsar light curves as observed in different energy ranges (starting from radio up to 
sub-TeV). This information will allow us to constrain relative location of different 
emission regions within the pulsar magneto sphere. Note that some theoreticians 
postulate the existence of a few emission regions (acceleration gaps) which might 
be responsible for the production of different parts of the pulsed spectrum. There 
is some evidence of the existence of the third peak in the pulsar's light curve, e.g. 
the case of Vela pulsa r [114 1. see also some features in the Crab light curve at sub- 
TeV energies Possible discovery of such other components by the 

CTA will put in trouble present general models for y-ray emission. 

The localization of the acceleration and radiation regions within the pulsar 
magnetospheres (thanks to the precise sub-TeV measurements) will make possible 
a more precise modelling of the e ± pair cascade processes and the e ± pair content 
in the inner pulsar magneto sphere. Note that the present models are not even able 
to give satisfactory predictions of the numbers of e pairs escaping through the 
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pulsar light cylinder into the PW Ne. T he calculations of e ± pair multiplicities in 
the inner pulsar magnetospheres lll8lll are at least an order of magnitude lower 
than their estimates based on the radiation output from their parent PWNe [|ll8ll . 
Sub-TeV y-ray observations of the pulsed emission are also expected to constrain 
a possible role of ions in the pulsar magnetospheres (by limiting the regions of 
electron acceleration) which are expected to carry the currents necessary to sup- 
port the pulsar acceleration and radiation mechanism. 

It will be also very interesting to study the stability of the sub-TeV pulsed 
component over long periods of time. Are the radiation processes responsible 
for the y-ray emission related in some way to the pulsar glitches or the giant 
radio pulses? The stability of sub-TeV emission and its possible relation to the 
variability observed in the radio will require investigation of the y-ray emission 
on short time scales. Therefore an instrument with an order of magnitude better 
sensitivity (like CTA) is very welcome. 

2.2. Millisecond pulsars 

Millisecond pulsars had not been widely expected to produce observable pulsed 
y-ray emission due to a much lower number of them observed in radio, their rel- 
atively low energy loss rates (much older objects, with weak surface magnetic 
field), preferable occurrence in binary systems, and theoretical predictions. There- 
fore, the discovery of a large number of MSPs by Fermi-LAT was rather a surprise. 



Before launching of the Fermi and AGILE telescopes 1120 ill pulsed y-rays were re- 
ported from only one MSP, PSR J02 18+4232. 

2.2.1. Discovery of the MSP population 

At first, 8 MSPs were discovered by using radio ephemeris [Q. Since MSPs 
are relatively close objects, many of them were expected to be found at high galac- 
tic latitudes. In fact, radio searches revealed several more MSPs in the Fermi-LAT 
error boxes of sources outside the Galactic disk. They show spectral properties 
typical for classical pulsars (i.e. the characteristic cut-off at a few GeV). The 
y-ray emission features of the MSPs resemble very much those observed in the 
case of CPs what suggests a common mechanism for y-ray production. In spite of 
significant general differences between these populations of pulsars (mentioned 
above), some other basic properties, e.g. the strength of the magnetic field close 
to the light cylinder, are very similar. This supports arguments that the emission 
region in MSPs is also located far away from the pulsar surface and the radiation 
process is similar. 



9 



Some of the discovered MSPs are in binary systems containing underluminos 
comp anions (the so c alled Black Widow pulsars) or typical low mass companions 



(e.g. 111381, 1193L I248I1. Many MSPs have also been discovered within Globular 
Clusters (GCs). Their cumulative emission is probably responsible for the ob- 
served y-ray emission from several GCs (as discussed in more detail in Sect. 4). 
In fact, in the case of one MSP within GC, PSR J 1824-2452 in M28, the pulsed 
y-ray emission has been already reported fl235J]. 

2.2.2. y-ray production models: expectations and observations 

In the past, MSPs were expected to be better sites for production of y-rays 
with energies detectable by Cherenkov telescopes than classical pulsars due to 
significa ntly weaker sur face magnetic fields and the resulting smaller y-ray ab- 



sorption 111221 13 131 1173] . The models discussed in those papers have predicted 
curvature y-ray spectra extending up to a few tens of GeV, and also a much weaker 
separate IC component extending up to ~1 TeV, produced by a comptonization of 
the thermal X-ray emission from the NS surface. The expectation of the curva- 
ture emission extending up to a few tens of GeV have not been confirmed by the 
Fermi-LAT observations. 

The calculations of the y-ray spectra from MSPs, performed in terms of other 



models, showed breaks at lower energies (e.g. [294]). This is due to the farther 
location of the y-ray production region from the surface of the neutron star. These 
results are more consistent with the present Fermi-LAT observations. Recently, 
calculations of the y-ray light curves obse rved from MSPs were made in terms of a 
few different models by Venter et al. [|296|| . It is concluded that most MSPs require 
a location of the emission region in the outer magneto sphere, what is consistent 
with the prediction of the outer gap and slot gap models. This conclusion is similar 
to that reached from the analysis of the emission models of the CPs, thus opening 
the possibility of a detection of sub-TeV emission also from the MSPs. 

2.2.3. Future CTA observations 

It is expected that the next generation Cherenkov telescopes will be able to 
establish possible sub-TeV y-ray emission from MSPs. Such second y-ray com- 
ponent in the pulsed spectrum, produced in the IC process, will be interesting 
to be investigated with the present and future Cherenkov telescopes. Does this 
emission extends with a simple power law after the peak at a few GeV (as seems 
to be observed already in the Crab pulsar), or is it a new separate component as 
expected in the earlier models for the MSP y-ray emission? Note that due to the 
large differences in the location of the light cylinder radius in CPs and MSPs, 
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the thermal radiation from the NS surface of the latter might be of much greater 
importance. In case of a positive detection, comparative studies of the sub-TeV 
y-ray emission features of CPs and MSPs (such as y-ray light curves, synchro- 
nization of the sub-TeV emission with emission in other energy ranges, cut-offs 
in the energy spectra), will provide stringent constraints on the location of high 
energy processes in the pulsar magnetospheres with different magnetic fields and 
rotation periods. 

2.3. Magnetars 

Magnetars (AXPs and SGRs) are rotating neutron stars with the surface mag- 
netic field about two orders of magnitude larger than that observed in classical 



pulsars [156]. On the other hand, their rotational periods are much longer being in 
the range 2-16 s. It is expected that about 20% of all neutron stars have been born 
as magnetars. Since their X-ray luminosity is clearly larger than the rotational 
energy loss rate, it has been proposed that these sources are powered by the decay 
of the super- strong magnetic fields in contrast to the MSPs and CPs powered by 
the rotational energy losses. Magnetars are relatively young objects, with the age 
comparable to the observed classical y-ray pulsars, i.e. 10 3 — 10 5 years. 

2.3.1. Expectations of pulsed y-ray emission from magnetars 

The processes occurring within the magnetospheres of magnetars are expected 
to show general similarities to those occurring in the classical pulsars in spite of 
large differences in the surface magnetic fields. In fact, the values of the magnetic 
field strength close to the light cylinder in magnetars are of the same order as those 
in CPs. Therefore, the radiation models developed for CPs (especially their outer 
atmosphere versions) have been applied also for magnetars. The calculations per- 



formed in terms of the outer gap model 111 341 131211 show that magnetars are likely 



sources of GeV y-rays detectable by the Fermi-LAT telescope. These theoretical 
expectations were recently confirmed by independent calculations 

However, recent observations of a few magnetars do not confirm these predic- 
tions. The uppe r lim its are clearly below predictions at least in the case of some 



magnetars (see [|262ill6ll . This requires either a re-consideration of the applicabil- 
ity of the outer gap model for the magnetar magnetospheres or a consideration of 
another scenarios for the magnetar phenomenon. A way to understand the prob- 



lem is the model of twisted magnetar atmospheres fllOU llOOtl . In this model e ± 
pairs are copiously produced already close to the surface of the neutron star as 
a result of magnetic e ± pair production by y-rays very close to their production 
place. These pairs limit the available potentials within the magnetosphere to ~10 9 
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V. In such case, only hard X-rays to soft y-rays (with very hard spectrum) can be 
produced by leptons which scatter re s onan tly thermal X-ray photons (in consis- 
tence with the observations ||202 , 169 , 160 1. 



Note however, that there seems to be a continuum transition of properties be- 
tween typical magnetars (B > 10 14 G) and CPs. Recently, GeV y-ray emission 
was discovered from a strongly magnetised, rotation powered pulsar (PSR Jl 1 19- 
6127) with the surface magnetic field 4.1 x 10 13 G, i.e. very close to the critical 
value B CL = 4.4 x 10 13 G when the electron cyclotron energy equals to its rest 



mass 1123 811 . In this case, quantum effects have to be taken into account when 
considering radiation processes close to the NS surface. This may significantly 
influ ence the possible development of cascades in the pulsar magnetosphere (see 



e.g. fl 1761 126411 . However, the y-ray emission from PSR Jl 119-6127 has features 



similar to those observed in CPs, again arguing for its origin in the outer magne- 
tosphere where the magnetic field is relatively weak. 

Halzen et al. fll V2h tried to estimate the fluxes of TeV y-rays and neutrinos by 
simple extrapolation from the hard-X-ray emission observed from magnetars. The 
y-ray fluxes predicted in these calculations are clearly above the present observa- 
tional limit reported by the VERITAS Collaboration for two magnetars [|l70ll . 
Note that the magnetic field strength close to the light cylinder in magnetars, clas- 
sical and millisecond pulsars is comparable. Therefore, it makes sense to expect 
pulsed y-ray signals also from magnetars. This pulsed emission, if originated in 
the outer magnetosphere, could extend up to sub-TeV energies as recently discov- 
ered in the case of the Crab pulsar. 

2.3.2. Can CTA detect y-ray emission from magnetars? 

Since the magnetic field strengths at the light cylinder in different types of ro- 
tating neutron stars are similar, it makes sense to search for pulsed sub-TeV emis- 
sion also from magnetars. The old arguments about the saturation of the electric 
field in the inner magnetosphere of the magnetar as a result of y-ray absorption in 
a strong magnetic field may not be so convincing, as it was previo usly expected, 
due to the formation of the positronium stage by the y-ray photons ll264ll . 

The rotational periods of magnetars were likely much shorter at birth than ob- 
served at present. As classical pulsars, they probably produce relativistic winds 
which terminate producing shock waves. Therefore, magnetars might also be 
surrounded by nebulae filled with relativistic particles as observed in the case of 
young pulsars. In fact, first extended TeV y-ray source has been recently discov- 
ered in the direction of the well known magnetar SGR 1806-20 [|256ll . A confir- 
mation of the association of this TeV emission with magnetar by CTA would be of 
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Figure 2: Comp ilation of the TeV y-ray measurements of the steady emission from the Crab Neb- 
ula (from [308], reproduced by permission of the authors). Note the localization of the maximum 
in the spectrum (cumulative MAGIC and Fermi data) and the uncertainties in its shape above a 
few TeV. 



great importance for comparative studies of the properties of the relativistic winds 
from magnetars and classical pulsars and for the constraints of the multiplicity 
of e ± pairs produced around NSs under very different conditions. Note that the 
nebulae around magnetars are expected to be more extended due to a larger en- 
ergy output at the early time after their formation. They should be characterised 
by weaker magnetic fields in comparison to the nebulae around classical pulsars. 
Therefore, they might resemble more the unidentified TeV sources which are eas- 
ier to discover in the TeV y-ray s than in the X-ray or radio observations. 



3. Pulsar Wind Nebulae 



The nebulae around energetic rotation powered pulsars have been suspect ed as 



sites of particle acceleration and y-ray emission since their discovery (e.g. [168 



250]. The high energy emission from these objects is expected to be produced by 



extremely energetic electrons interacting with the magnetic field and background 
radiation within the nebulae. Up to now, several such nebulae have been reported 
to emit TeV y-ray s. A variety of their spectral and morphological features have 
been discovered. At present, however, they can not be investigated in enough 
detail due to the limited sensitivity, angular and energy resolution of the present 
telescopes. 
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3.1. Observations of the Crab Nebula 

The Crab Nebula was the first object of thi s type which was established as a 
TeV y-ray source about twenty years ago 00311 . Its multiwavelength spectrum is 
composed of two broad bumps wildly interpreted as the effect of the synchrotron 
and IC radiation of electrons with energies up to a few 10 15 eV (e.g. jl43L \72\ . 
The spectrum showed an amazing stability within a few percents in most energy 
ranges. Therefore, the Crab Nebula has been considered as a standard source 
in the y-ray astronomy 022011 . However, some variability at the end of the syn- 
chrotron spectrum (in the y-ray range) on a time scale of a few mon ths was men- 
tioned based on the EGRET and COMPTEL observations rf226L fl44ll . Apart from 
these observations, the broadband spectrum of the Crab Nebula seemed to be well 
known and steady except the energy range above a few TeV where some differ- 
ences might be noted resulting probably from the lack of statistics, see observa- 
tions by the HEGRA H, HESS H, and MAGIC telescopes Hi]. The Crab 
Nebula spectrum reported by Milagro in the energy range 1-100 TeV is generally 
consistent with the previous HERGA measurement Q23D . For a collection of these 
high energy observations see Fig. 2. Note that recently a long time scale, low 
level variability of the hard X-ray emission from the Crab nebula was reported 
(e.g. H. 



It came as a big surprise when at first the AGILE team [|283|l and next the 
Fermi-LAT team [|2l|l reported the discovery of two large flares (flux increase 
by a factor of 4-6) lasting for several days at energies above 100 MeV from the 
Crab Nebula (see Fig. 3). More detailed analysis of the September 2010 flare 
shows separate spikes on a half-day time scale 117 611 . Even more drastic flare 
from the Crab Nebu la was observed in April 2011 (see Fig. 4, the details re- 
ported in Ill20l 127211 ). During this flare the emission increased by a factor of 
~30 on an hour time scale. The averaged spectrum of this flare has been well 
described by a pow er law with the spectral index 1 .6 and the exponential break 
at ~580 MeV O- 

Moreover, no accompanying variability in the other energy 
ranges have been observed, except the report on a simultaneous enhancement of 
the signal at ~1 TeV by the ARGO-YBJ Collaboration ||5fl], not confirmed by 
MAGIC VERITAS (Si or MILAGRO jivfo . 



3.1.1. Interpretion of the variable y-ray emission from the Crab Nebula 

These observations pose serious problems for the popular mechanism of parti- 
cle acceleration at the pulsar wind shock since they seem to violate the constraints 
on the maximum energy of accelerated particles, obtained by assuming diffusive 
shock acceleration or any mechanism connected to acceleration in a magnetic field 
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on a timescale longer than the gyration timescale. Note that this limitation can be 
overcome in the case of the magnetosonic wave interaction with the termination 
shock [2131]. Moreover, the emission region is expected to move relativistically 
towards the observer (as considered in [93 J) and/or the acceleration of particles 
can occur in other acceleration mechanism such as e.g., reconnection of the mag- 
netic field in the pulsar wind [93, 293ll or in the large scale electric field generated 
by the rotating pulsar ll2lll . In another model, a variable synchrotron emission is 
proposed to be produced due to stochastic variations of the magnetic field in the 
reconnection region (112311 . The acceleration site (the place of the flare origin) is 
also enigmatic. It has to be very compact, i.e. much smaller than the travel time 
of light across the shock region. It ha s to b e identified with a small region within 
the pulsar wind, e.g. inner knot? (see 

This surprising phenomenon will certainly reach much attention of the CTA, 
HAWC and other Cherenkov telescope projects. The Crab Nebula monitoring at 
multi-TeV energies by the water Cherenkov detectors seems to be of great im- 
portance since the flaring synchrotron component is expected to be produced by 
the ~PeV leptons which might eventually also produce flaring Inverse Compton 
emission at multi-TeV energies. This can only occur when some specific condi- 
tions are fulfilled, i.e. leptons can cool on the ICS on the time scale of the flare. 
This process can become efficient in the case of a significant Doppler boosting of 
the acceleration region of leptons V93n . On the other hand, a possible discovery 
of variable multi-TeV emission from the Crab Nebula with the next generation 
instruments would allow us to estimate the maximum energy of the accelerated 
leptons and thus constrain the emission region of the flaring component. Such 
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Figure 4: The spectra of flaring y-ray emission (see Fig. 3) from the Crab Nebula during the de- 
velopment of the flare in 201 1 (from lll20lD . Note the behaviour of variable synchrotron spectrum 
observed at different moments of the flare and the maximum of the spectrum extending up to ~1 
GeV. (Reproduced by permission of the AAS) 



variability would be in contradiction with the model U2M- 



3.1.2. Observations of other PWNe and their interpretation 

Observations of other PWNe at high energies show a variety of morphological 
and spectral features. In general, it is not so surprising since pulsars responsible 
for the PWNe are expected to be born with different initial parameters (period, 
surface magnetic field, velocity), sometimes in an inhomogeneous dense medium 
still present around them due to relatively short evolution time of their progenitor 
stars. The parameters of a newly born pulsar determine the total energy content 
of the future nebula, the magnetic field inside it and also (likely) the maximum 
energy of the accelerated particles. On the other hand, the surrounding medium 
determines the size of the nebula, its morphology, and maybe also the dominant 
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type of y-ray radiation mechanisms (e.g. in some cases bremsstrahlung from elec- 
trons and pion decay from hadronic interactions). Recent observations, mainly in 
the radio, X-ray and TeV energy ranges, reveal this variety of features. However, 
the sensitivity and angular resolution of the present TeV y-ray telescopes do not 
allow us to investigate of the details of the spectral and morphological features in a 
large population of objects (only in exceptional objects energy dependent morpho- 
logical studies have been performed). Therefore, robust conclusions about their 
physics can not be reached at present. Below we will review the most interesting 
features observed in some specific objects. 

A very extended TeV nebula (~ 0.5° x 1.5°), Vela X, has been discovered in 
the vicinity of the well known Vela pulsar. It is clearly shifted with respect to 
the pulsar location. The spectrum of the Vela X nebula is flat peaking at a few 
TeVs, atypical in comparison to other PWNe Q39[|. On the other hand, recently 
discovered GeV emission from the direction of Vela X 11171 l236ll is much more 
extended, showing a spectrum with completely different features (see Fig. 5). The 
emission from a more compact TeV nebula (cocoon) was modelled assuming its 
hadronic Ill86ll or leptonic origin [|203ll . The complicated structure of the Vela 
nebula in different energy ranges was explained as due t o the two populations of 
electrons injected by the pulsar at present and in the past 

Other PWNe, observed in the GeV and TeV energies, also show spectra with 
clear discontinuities suggesting an existence of two populations of particles, e.g. 
HESS J1640-465 |40, l27lll . In the case of young objects this might be a result 
of the acceleration of particles within the PWNa itself and also on the supernova 
shell, sometimes visible in the so called composite PWNe. The two component 
spectra are not characteristic for all nebulae. For example, the nebula MSH 15-52 
shows the GeV- TeV spectrum nicely described by the IC emission from a single 
population of electrons, and the radio to X-ray spec trum - by the synchrotron 
process from this same population of electrons [ 261 , [3~3ll . As in the case of the 
Crab Nebula, the GeV spectra of other PWNe seem to be clearly harder than the 
TeV spectra (e.g. MSH 15-52 or PSR J1023-5746/HESS J1023-575 Q- 
They are produced by electrons with lower energies characterised by long cooling 
time scale within the nebula. Their investigation should provide information about 
the injection rate of electrons in the early age of the nebula. Therefore, precise 
measurements of the relative emission from nebulae in the GeV and TeV energy 
ranges should allow us to estimate the total power in relativistic particles and so 
to constrain the initial parameters of the parent pulsars (see e.g. Ill46ll . 

The morphological features of PWNe in the TeV and X-ray energies make 
possible to distinguish two groups which differ in the age. Young PWNe show 



17 



~ io- 10 



Halo 




Cocoon 



s 

o 

M 

10 " U 

Ed 
■O 
\ 
Z 
T3 

H 10 -ia 



5 

log Energy (eV) 



i 1 1 1 1 1 


1 1 I 1 I 1 1 1 I 


1 1 1 1 1 1 i 




/"""V ASCA 


HESS : 


> 


1 \ , , 





10 



15 



log Energy (eV) 

Figure 5: Compilation of spectral measurements of different regions (upper panel - Halo and lower 
panel - Cocoon) in Vela-X pulsar wind nebula collected from different instruments (for details 



see B17I1 . reproduced by permission of the AAS). Note large differences in the spectral shape 
between the Cocoon (produced by fresh particles from the Vela pulsar) and the Halo (produced by 
old, lower energy particles). 
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compact (X-ray and TeV) nebulae with the pulsar more or less in the middle. On 
the other hand, older PWNe are extended with the pulsars not located centrally. 
They are clearly more compact in X-rays than in TeV y-rays. This feature has been 
expected since synchrotron X-ray emission is produced by electrons with larger 
energies in the regions closer to the pulsar where the magnetic field is stronger. 
On the other hand, the TeV emission comes from the electron IC scattering on 
the background radiation in the vicinity of the PWNa (such as the Cosmic Mi- 
crowave Background (CMB) and the infrared-optical background in the Galactic 
disk). These general features are consistent with the recent measurements of the 
structure of the PWNa, HESS J1825-137, showing a softening of the TeV spec- 
trum with the increasing distance from the pulsar [14111 . 

The complex structure of the PWNe may be a consequence of its interaction 
with the asymmetric reverse shock of the supernova remnant [106] . This can result 
in a displacement of the nebula with respect to the pulsar as frequently observed 
in TeV energies. Also the high velocity pulsars can result in a formation of asym- 
metric PWNe. What influences the morphology of a specific nebula requires an 
investigation of its detailed structure in the y-rays with a sensitive instrument and 
also measurements of the pulsar velocity and the PWNa surroundings. For exam- 
ple, in the case of the Vela PWNa, the asymmetric structure of the nebula with 
respect to the location of the pulsar is likely due to the interaction with the reverse 
shock since the pulsar moves in a completely different direction than the largest 
extent of the nebula [jl50 1. 

The structure of the PWNa may not only be determined by its background soft 
radiation content and the surrounding medium (including supernova remnant) but 
also it should depend on the geometr ical structure of particle injection. It is wel l 
known from X-ray observations (e.g. II 1 8011 and MHD simulations (e.g. II109lI301I1 
that the acceleration process of electrons is likely to occur highly anisotropically 
(e.g. the presence of the equatorial pulsar wind shock, jets, etc.). Probably we 
are in touch with such a case in the nebula MSH 15-52 around the pulsar PSR 
1509-58, where the TeV nebula is aligned with the direction of the X-ray jet [|33ll . 

Possible presence of dense clouds in the PWNa vicinity might result in addi- 
tional components in the GeV-TeV spectrum due to the interaction of electrons 
with matter (bremmstrahlung process) or possible interaction of hadrons (pion 
production). Such interacting PWNe might be good indicators of the accelera- 
tion of hadrons by the PWNe 117911 . Detailed spectral and morphological studies 
with CTA and in other energy ranges should answer the question whether such 
objects exist at all, what is the relative contribution to the high energy range from 
the PWNa, supernova shell, and particles which escaped from the PWNa and the 
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supernova shell. 



3.1.3. How CTA can help to understand PWNe? 

Even in the case of the best studied objects, such as the Crab Nebula, CTA 
with its angular and spectral resolution can answer some new enigmas. CTA 
is expected to contribute to the problem of the extent of the TeV source in the 
Crab Nebula. Is it really a point like source (with a dimension corresponding to 
the pulsar wind shock) or is it comparable to the extent of the whole nebula (2- 
3 arc min) or maybe even larger? Is it possible to constrain energy dependent 
dimension of the Crab Nebula emission? Very important question is related to the 
variability of the Crab Nebula multiwavelenght spectrum. Can the whole spectrum 
above several GeV (Inverse Compton origin) vary on a long time scale as recently 
discovered in the X-rays? It would be very important to establish the spectral 
features at the highest energy end of the Crab Nebula spectrum (i.e. at ~10 TeV). 
Is there any evidence of a contribution from an other radiation process (not IC) as 
postulated by some hadronic models (e.g. 11971. l64l 19111 ? This part of the spectrum 
could also change on short time scale in accordance to the variable emission above 
~100 MeV since both are expected to be produced by electrons with energies of 
a few hundred TeV. A determination of the end of IC y-ray spectrum from the 
Crab Nebula by CTA would also allow us to constrain the magnetic field strength 
in the inner nebula since the location of the cut-off in the TeV region should also 
depend on the relative efficiency of the lepton cooling on the synchrotron and IC 
processes. 

The very large angular dimensions of some nearby PWNe (even of the order 
of a few degrees) will require an instrument with a very large field of view. CTA, 
with a field of view even larger than that of the present HESS field of view and an 
order of magnitude better sensitivity, will make possible detailed studies of fine 
structures within PWNe. Note that old close nebulae are expected to contain pul- 
sars which produce bow shocks due to their fast motion through the interstellar 
medium. On such bow shocks electrons might be accelerated up to the TeV ener- 
gies. A few of such bow shocks have been observed up to now in the X-rays (e.g. 
around PSR J1747-2958, lfl65ll ) or around Geminga, e.g. ll239ll ). The low level 
TeV emission from the bow shock, together with the large scale emission from 
the PWNa, produced by older electrons injected in the vicinity of this same pul- 
sar, would enable us to make comparative studies of the two y-ray objects created 
by the same pulsar. Detection of such PWNe with CTA would provide stringent 
constraints on evolution models of the PWNa and on the pulsar wind content. 

Detailed investigation of energy dependent morphology of such composite 
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PWNe with the future CTA will certainly allow to conclude on the location of 
the dominant particle acceleration site within such nebulae. In the case of such 
nebulae particles can be accelerated not only in the vicinity of the pulsar (pulsar 
wind shock, bow shock) but also by the supernova shocks (forward and reverse) 
appearing as a result of their interaction with the inhomogeneous surrounding 
medium. 

It will be very interesting to study the spectral and morphological features 
of the Vela nebula from sub-TeV up to multi-TeV energies in order to find out 
whether the y-ray emission has actually either a two component origin or it is due 
to a single component which changes continuously its properties with the distance 
from the center of the nebula. Only a telescope with the sensitivity and angular 
resolution of CTA can put new light on this problem. 

What factors are responsible for the variety of the spectral features observed 
in PWNe remains a mystery. This problem can be answered with accurate ob- 
servations with CTA provided that it has good sensitivity already at ~20-30 GeV 
and also above ~30 TeV. The good sensitivity at large energies should allow us to 
constrain the maximum energies of accelerated particles. The information on the 
end of the electrons spectrum and on the end of synchrotron spectrum (expected in 
the hard X-ray - soft y-ray energy range) will allow us to independently constrain 
the magnetic field strength within the PWNe. 

The radiation efficiencies of the PWNe are estimated in the range of a few up 
to ~ 30%. They constrain the so-called sigma parameter describing the ratio of 
energies injected into the nebula in the form of the magnetic field and relativistic 
particles. The knowledge of this value and its evolution with the pulsar age is 
the basic factor limiting the multiplicity of pairs injected from the inner pulsar 
magnetosphere and the acceleration processes occurring in the pulsar wind region 
and/or the pulsar wind shock. Therefore, reliable measurement of the radiation 
efficiency in the large population of PWNe with CTA will be very important for 
constraining the pulsar models determining those of the PWNe. 

CTA should also have a possibility to resolve fine structures (e.g. peaks?) in 
the TeV y-ray spectrum which might be due to direct comptonization of the soft 
radiation by leptons in the relativistic pulsar wind [110811 . Their discovery with 
CTA would indicate to the region of the relativistic wind formation and help to 
determine its Lorentz factor. We may wonder whether this Lorentz factor scale 
in a simple way with the basic parameters of pulsars and how it depends on the 
evolutionary stage of the PWNa. 
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4. Globular Clusters 



Globular Clusters (GCs) are concentrations of 10 5 — 10 6 old stars (with the age 
> 10 Gyrs) contained within a spherical volume of a few pc. About 150 clusters, 
discovered up to now, create a spherical halo around the Galaxy with typical radius 
of ~10kpc Ill77n . GCs also contain plenty of compact objects, millisecond pulsars 
(MSPs), Cataclysmic Variables (CVs) and Low Mass X-ra y Binaries (LMXBs). 
Up to now, about ~140 MSPs have been discovered in GCs lfl24l l 16311 . and thou- 
sands of X -ray sources are observed. They are expected to belong to the class of 
CVs (e.g. 1129 811 ). It is argued that the numbers of these compact sources (MSPs, 
CVs, LMXBs) in a specific GC correlate with the so called enco unter rate which 
is the combination of the core density and the core radius of GC Il245l.l246l. il 811. 



4.1. Gamma-ray observations 

Several Globular Clusters have been recently established as sources of GeV 
y-rays in the observations with the Fermi-LAT telescope [jl, 199 1. They are char- 
acterised by differential spectral indices in the range 0.7 4- 1.4 and exponential 
cut-offs at 1.0 4-2.6 GeV. The y-ray power of various GCs differ by an order of 
magnitude. Energy spectra of GCs show features very similar to those recently 
observed in the population of the MSPs [18]. However these features may not be 
common for all GCs since some of the recently discovered objects show spectra 
extending above ~10 GeV, without clear evidence of cut-offs |275j]. Therefore, 
the origin of y-rays may not be only related to processes occurring within the inner 
magneto spheres of MSPs. 

The TeV emission has been searched from GCs a lread y with the Whipple tele- 
scope (see the upper limits for M13 and M15 IU7U I206I1 ). More stringent upper 
limits, at the level of present theoretical predictions, wer e ob tained by modern 
Cherenkov telescopes for M13 H, M5, Ml 3 and M 15 ll2~18ll and Tuc 47 & 
Only one TeV y-ray source was discovered close to the GC, Ter 5 112411 . Ter 5 con- 
tains the largest number of the MSPs (33 objects). The observed y-ray flux (1.5% 
of the Crab Units) is at the level expected from the models but the source is clearly 
shifted from the position of Ter 5 . On the other hand, the position of the TeV y-ray 
source is consistent with the location of the diffuse nonthermal X-ray source (115811 
and the radio source fll37ll . However, why is the source shifted from the position 
of Ter 5 is at present unknown. The energy spectra observed from the direction of 
Ter 5 in the GeV and TeV energies can not be connected smoothly (see Fig. 6). 
They seem to constitute two independent components. This suggests their origin 
by different radiation mechanisms in different locations, e.g. GeV emission from 
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the GC itself and the TeV emission from a completely different type of source, 
maybe even not directly related to the GC itself. 

The y-ray emission from GCs shows some interesting correlations with the GC 
content. For example, Abdo et al. [[6] and Hui et al. [il87ll find a clear correlation 
of the observed GeV y-ray luminosity with the stellar encounter rate which deter- 
mines the number of MSPs and Cataclysmic Variables within the specific GC. It is 
believed that the encounter rate determines the formation rate of compact binaries 
which are responsible for the origin of the MSPs inside GCs. The encounter rate 
can be estimated from the observed density of stars in the core of GC and the core 



radius [I166I1 . Moreover, Hui et al. Hi 8711 report the positive correlation of the y-ray 
luminosity with the metallicity [Fe/H] of the GC. They also find a tendency of the 
correlation of the y-ray flux with the energy density of soft photon field at the GC 
location. 

The GCs identified with the Fermi-LAT y-ray sources differ in their GeV y- 



ray power by an order of magnitude [18]. Based on the observed y-ray luminosity 



of a specific GC, the average energy loss rate of MSP (adopted for MSPs ob- 
served in Tuc 47 [@], and the average spin-down to y-ray luminosity conversion 
efficiency [@, [Tsb allows one to estimate the number of MSPs in specific GC. 
The populations of MSPs in specific GCs derived by Abdo et al. [| 1 IN are in good 
agreement with the estimates of the number of MSPs obtained with the help of 



other methods (see the estimates for Ter 5 in II1641 [199]) 



4.2. Theoretical predictions and confrontation with observations 

Since globular clusters contain many MSPs, it was expected that their cumu- 
lative GeV-TeV y-ray emission might be above the sensitivity of the y-ray tele- 
scopes. A few estimates of the y-ray fluxes from MSPs in globular clusters were 
published already before the launching the Fermi Observatory. Based on the pair 
starved polar cap model Harding et al. predicted the y-ray fluxes from Tuc 47 
near the EGRET upper limits. Venter & de Jager (129511 . applying the modern ver- 
sion of the polar cap model (general relativistic effects included), also predicted 
y-ray spectra well above the sensitivity of the Fermz'-LAT telescope. However, 
there are at present some doubts concerning the validity of the polar cap model as 
a likely scenario for the y-ray emission from the pulsar inner magnetospheres (see 
Section 2). 

The presence of a very dense radiation field within the clusters, produced 
by normal stars (2-3 orders of magnitude larger than the energy density of the 
CMB), suggests that relativistic particles injected by MSPs (and other compact 
sources) should efficiently produce TeV y-rays in the Inverse Compton Scattering 
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Figure 6: The 7-ray spectra observed from the Globular Cluster Ter 5 in the GeV energy range 
(left figure from IU8I1 . reproduced with permission of A&A) and in the TeV energy range (right 
figure from 024I1 . reproduced by permission of the A&A). Note that the TeV component has likely 
different origin than the GeV component since the spectrum index, estimated from the two last 
points of the Fermi spectrum, is much steeper than the spectral index observed in the TeV energy 
range. 



process. Such general scenario was considered for the first time by Bednarek & 
Sitarek Ir9~8ll . In this model, it is assumed that leptons are accelerated either in the 
millisecond pulsar winds or shocks produced by such winds (as a result of their 
collisions with the stellar winds or collisions of the two pulsar winds between 
themselves) or they are injected directly from the inner pulsar magnetospheres. 
These energetic leptons diffuse through the volume of the GC and interact at the 
same time with its soft radiation, i.e with the optical photons from the stars and 
the CMB. Due to the lack of detailed knowledge about the spectral features of 
leptons injected by the MSPs, a few different models for electron injection are 
considered with the limitations expected from a comparison with the classical 
pulsar population. The basic properties of the y-ray spectra calculated in the case 
of mono-energetic injection of leptons into the soft radiation field dominated by 
soft photons from the stellar population within the cluster and the CMB are shown 
in Q98U . In terms of such a model, the y-ray spectra were calculated for selected 
globular clusters assuming different spectra of injected leptons [98] . A compar- 
ison of the predicted spectra with the observations of specific GCs allowed the 
authors to put the upper limits on the product of the number of MPSs within GC 
and their lepton injection rate. 

Similar models were also considered in [|296 , 1 33h . In the first paper it is as- 
sumed that leptons are injected from the inner magnetospheres of the millisecond 
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pulsars through the light cylinder in the pulsar polar cap model [129511 . In such a 
case, the injection rate of leptons is directly linked to the pulsed y-ray emission 
from the MSPs observed in the GeV energies. The TeV emission expected in 
this work is quite similar to that predicted by Bednarek & Sitarek. On the other 
hand, Cheng et al. II 13311 consider a production of y-rays in the IC scattering on 
the galactic infrared and optical background by leptons in the pulsar winds. They 
predict an appearance of a relatively narrow second bump in the y-ray spectrum 
from GCs at energies below ~100 GeV. However, the abrupt cut-off in this spec- 
trum above 100 GeV seems to be in contradiction with the observations of the 
TeV emission from Ter 5, provided that this emission comes from the GC itself. 

It has also been suggested that a significant part of the TeV emission ob- 
served from GCs may also come from other sources within the GC. For example, 
Bednarek Il90ll notes that a typical GC should contain several thousands of WDs 
(with a significant fraction with the surface magnetic field of the order of 10 8-9 
G) which might be able to accelerate electrons to TeV energies (for more details 
see Section 9.3). On the other hand, Domainko lll51[l and Abramowski et al. 112411 
propose that a recent type la supernova explosion or a short Gamma-Ray Burst 
can accelerate hadrons which could produce y-rays in the interactions with matter 
within or around the GC Ter 5. 



4.3. Future studies of GCs with CTA 

Sensitive observations of GCs in the sub-TeV energies (below ~100 GeV) will 
be very important for a determination of a relation between the observed spectra 
of the Fermi-LAT sources and the predicted TeV emission. Do GCs really emit 
TeV y-rays produced by the leptons accelerated by the relativistic winds of the 
MSP population? 

The next generation of Cherenkov telescopes should clarify the origin of the 
TeV source in the vicinity of Ter 5. Are there any morphological signatures of the 
relation of this source to Ter 5? This problem could be answered by the energy 
dependent morphological studies of the HESS source towards Ter 5. For example, 
the closer geometrical connection between the HESS TeV source and the center of 
Ter 5 in the multi-TeV energies than in the sub-TeV energies will suggest that the 
relativistic particles originate within this GC. Detailed investigations of spectral 
and morphological features might also allow us to assign this source to the already 
known classes of the TeV sources (SNRs, PWNe, binary systems)? 

Detection of sub-TeV pulsed emission from specific MSPs within the GCs 
would make possible to estimate relative contributions to observed GeV emis- 
sion, with respect to the total, of separate MSPs and that from the more extended 



25 



component produced by relativistic leptons escaping from their magnetospheres 
or perhaps also from magnetized WDs. This will allow to put additional con- 
straints on the parameters of the supposed relativistic winds from the MSPs since 
they propagate in a well defined background radiation within the GCs. 

Another problem to be answered is related to the types of sources are poten- 
tially able to produce y-rays. Are the y-ray spectra of GCs similar to each other 
or do they differ significantly forming separate classes? At present it can not be 
excluded that the contribution to the observed y-ray emission from specific GCs is 
dominated by different objects expected to produce y-rays such as: MSPs, mag- 
netized WDs, Low Mass X-ray Binaries (LMXBs), type la supernova remnants 
or novae. There is also some evidence that the soft radiation content of several 



GCs can vary significantly [188]. It will be interesting to find out which of the 
two general types of GCs (centrally collapsed and extended) correlate in any way 
with their possible TeV emission. Only CTA with an order of magnitude better 
sensitivity can clarify these problems. 



5. Pulsars within binary systems 

For a long time X-ray binary systems have been suspected as sites of high 
energy processes which can lead to a production of y-rays. However, the first 
claims of the TeV-PeV y-ray detection from these sources were not confirmed by 
the more reliable observations with the Whipple Cherenkov telescope [302]]. In 
the GeV energies, a few X-ray binary s ystem s were detected in the large error 
boxes of the EGRET telescope (LS 5039 £3, Cyg X-3 (22! , LSI 61 303 lEil . 
Cen X-3 Q299[|). These associations were not conclusive since the y-ray signals 
did not show any evidence of the expected modulation with the orbital period of 
the binary systems. 

5.1. Recent observations in the GeV-TeV energy range 

An observational breakthrough has come with the next generation of the Cherenkov 
telescopes, which were able to discover clear TeV y-ra y sig nals from three massive 
binary systems, LS2883/PSR1259-63 H, LS5039 Q, and LSI 61 303 H. 
The signal in the first system was detected close to the periastron passage of 
the compact object which orbit around the massive star has well known param- 
eters H S- The observations of LS2883/PSR1259-63 with the Fermi-LAT 
telescope show quite unexpected features such as the brightening of GeV emis- 
sion about a month after the periastron passage [127611 . This GeV emission seems 
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Figure 7: The observations of the binary system PSR1259-63/SS2883 in the X-rays, GeV and 
TeV y-rays close to the periastron passage (from N276L reproduced by permission of the AAS). 
Note that the GeV emission (with the maximum at about a month after periastron) seems not to be 
correlated with the X-ray or the TeV emission. 



not to be directly correlated with that in TeV since its maximum appears at lat- 
ter period (see Fig. 7 for a collection of the observations close to the periastron 
of X-rays, GeV and TeV y-rays). LS5039 and LSI 61 303 showed clear modu- 
lation of the TeV signal with the periods of the binaries. The maximum emis- 
sion corresponded to the location of the compact objects in front of the massive 
stars 11541 14211 . The y-ray spectrum of LSI 61 303 is well described by a single 
power law without evidence of a spectral cut-off plQ. On the other hand, the 
spectrum of LS5039 is described by a single power law at the low state but by the 
power law with an exponential cut-off at the high state [42] . A clear correlation 
of the TeV and X-ray emission from the TeV y-ray binaries have been reported, 
thus supporting their production by the same populati on of electrons, e.g. LSI 61 
303 [66], LS 5039 Ell, and also HESS J0632+057 iDjLlh- These binaries have 
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Figure 8: On the left: Energy spectrum of LS 1 61 303 (Abdo et al. 2009f, reproduced by per- 
mission of the AAS). Note the clear deficiency of the y-ray emission at energies between ~10- 
100 GeV. On the right: Light curves of LS 1 61 303 reported by different Cherenkov telescopes 
(from [27], reproduced by permission of the AAS). There is no detection of the TeV emission at 
phases -0.6-0.8 by VERITAS in 2008-2010. 



been also detected in the GeV y-rays showing clear anticorrelation with the TeV 
light curve Isllsl]- It is know at present that the TeV emission, at least from LSI 61 
303, shows a long time scale variability durin g w hich the light curve changes sig- 



nificantly its shape (see Fig. 8 on the right H27U190I1 ). A comparison of the spectra 
measured from LS5039 and LSI 61 303 in the GeV and TeV energy ranges shows 
two distinct components, the first extending above and the second below ~10-100 
GeV (see Fig. 8 on the left lH,^. Such a feature was not expected by any theoret- 
ical model. Recently, two additional massive binary systems, HESS J0632+057 
detected in the TeV energies HQ] and 1FGL J101 8.6-5856 detected in the GeV 
energies [28], show modulation of the y-ray signal with the periods of their binary 
systems. 

5.2. Present understanding of the y-ray emission from binary systems 

The X-ray emission observed from the TeV binaries is atypically weak with 
respect to their y-ray emission, in contrast to the observations of other accreting X- 
ray binaries. Moreover, their radio morphologies suggest that the compact objects 



in these binaries are rotation powered pulsars [115211 . Therefore, the most popular 
scenario, which tries to explain the emission features from these objects, assumes 
an existence of a young energetic pulsar as already observed in the case of the 
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binary system LS2883/PSR1 259-63. The interaction of the pulsar wind with the 
strong wind from the massive star is responsible in this scenario for the formation 
of the shock wave within the binary system. Particles can be accelerated at this 
shock to at least TeV energies. Alternatively, they can be injected from the inner 
pulsar magnetosphere or accelerated in the pulsar wind region as expected in the 
case of isolated pulsars (see Sect. 2 and 3). 

Recently, a magnetar-like X-ray burst was observed from the direction of LS 
1 61 303, sugge sting that the pulsar can be characterised by a very strong surface 
magnetic field [289]. In fact, it was shown that particles c ould be accelerated to 
TeV energies also in such accreting magnetar scenario [1871 128911 . 



5.2.1. Leptonic models 

In the most popular scenario, y-rays in TeV binaries are produced by relativis- 
tic electrons which comptonize the thermal radiation from the massive compan- 
ion (e.g. 12151]). In fact, the optical depths within the massive binaries are clearly 
above unity for electrons with energies above a few tens of GeV Th us, h i gh e n- 
ergy y-rays are expected to be produced with large efficiency (e.g. [|247l . |224J]). 
However, in such a case their absorption becomes also important. The opti- 
cal dept hs w ere calculated for the massive stars in LS 5039 and LSI 61 303 by 
e.g. 11831. 1 14911 15311 . Due to large values of the y-ray optical depths, it is expected 
that the y-ray spectra from these binaries are formed in the IC e ± pair cascades 
propagating in the anisotropic radiation field of the massive star as already consid- 
ered some time ago by Bednarek Q80Q . The calculations of the spectra produced 
in the IC e pair cascade model, which es cape fro m these binaries, have been 



performed in several papers, see e.g. [84, |85l |43j [268 ] . These models predict clear 
anticorrelation of the GeV and TeV emission (e.g. Q83D) . This general feature has 
been later confirmed by the GeV-TeV observations of LS5039 and LSI 61 303. 
However, the basic conditions for the cascade process within the binary system are 
not well known. For example, the energy losses of the primary electrons and the 
secondary e ± pairs i n the casca de on the synchrotron radiation can play an impor- 
tant role (see H EH Hi • 

It was shown in these studies that, if the magnetic 
field in the cascade region is above ~ 1 G, a significant energy is transferred from 
leptons to synchrotron radiation. As a result, IC e ± pair cascade starts to be inef- 
ficient especially at TeV energies when the IC cross section in the Klein-Nishina 
regime drops significantly. By comparing the energy loss rates of electrons on the 
synchrotron and the IC processes, we can estimate the limit on the stellar mag- 
netic field strength below which the IC losses start to dominate. When the ICS 
occurs in the Thomson regime, it is simply given by the condition B < 407 4 2 G, 
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where T = lO 4 ^ K is the surface temperature of the companion star [18011 . For 
typical surface temperatures of stars in these binaries, of the order of 3 x 10 4 K, 
this critical value of the magnetic field is equal to ~ 400 G. Note however, that 
the energy density of magnetic field around the star drops in the interesting range 
of distances as t/s x R 1 4 , whereas the energy density of stellar radiation drops 
only according to Ub x R 2 . Therefore, we conclude that the region, where the 
synchrotron losses of electrons dominate over their IC losses, should be confined 
to that relatively close to the companion star. Provided that the cascade process 
does not occur deep in the Klein-Nishina regime, the synchrotron energy losses do 
not affect it strongly when the compact object is close to the apastron (the phase 
of the maximum of TeV emission), i.e. at the distance of 4 — 15 stellar radii. At 
these distances from the companion star, the magnetic fiel ds is like ly to be of the 



order of ~ 1 G, in accordance with the calculations in [80, 1113L 119411 



The magnetic field of a massive star can also change the directions of the 
primary and the secondary leptons resulting in different cascade scenarios. The 
following cases can be distinguished: (1) linear cascade in which particles (y-rays 
and leptons) follow the direction of the primary electron (e.g. tl26{]; (2) com- 
pletely locally isotropized cascade in which the secondary leptons are isotropized 
by the local magnetic field in their place of birth [|84[ |85[ 14311 ; and the most general 
case (3) magnetically driven cascade in which the secon dary leptons are driven by 
the local magnetic field with a certain structure fl266ll . Each of the above case 
predicts different angular distribution of the escaping y-ray photons on the sky, 
thus resulting in different light curves. Note that all these cascade scenarios can 
in principle be realized in different parts of the same binary system. 

Another important factor affecting on the escaping y-rays may be a variable, 
asymmetric, and clumpy wind from the massive star. A change of the stellar wind 
pressure should result in the change of the location of the shock within the binary 
system. Thus, the conditions for the acceleration of electrons and the development 
of the IC e ± pair cascade may change significantly, influencing the y-ray light 
curve. An aspherical wind, with its slow and dense equatorial part but fast and low 
density polar part, may result in a drastic change of the location of the shock with 
respect to the companion star at different phases of the compact object. This may 
strongly affect the conditions for the scattering o f the soft radia tion by electrons 
and the development of the IC e ± pair cascade ri26vL I269L fTToh . Note that also 
the pulsar wind is expected to be aspherical (e.g. QIOVD). Therefore, the geometry 
of the y-ray production in the IC e ± pair cascades initiated by electrons in the 
vicinity of the massive star may be very complicated. 

The winds of the massive stars in TeV binaries are expected to be clumpy. 
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Such dense clumps provide the medium for the propagation of relativistic parti- 
cles with very different properties (e.g. magnetic field strength, density of mat- 
ter). T his m ight result in a more complicated shape of the observed y-ray spec- 
trum r ibldl . Moreover, the y-ray light curve of the binary system can be also 
influenced by the relativistic flow of the plasma along the cometary tail of the 
pulsar on its orbit around the companion star [|l54ll . This effect is not expected 
to be very strong due to a rather slow pulsar wind behind the shock (velocities 
of the order of 0.3 — 0.5 of the velocity of light). The complicated shapes of the 
GeV-TeV spectra observed from LS 5039 and LSI 61 303 suggest the existence 
of two populations of electrons or two different emission mechanisms. In fact, 
an essential part of the GeV component could come from the inner pulsar mag- 
netosphere. On the contrary, the TeV component might be produced by electrons 
accelerated in the pulsar wind region 1128711 . In another scenario, the two popula- 
tions of relativistic electrons can be produced within the binary system as a result 
of the acceleration of particles on the two shock structure separated by a contact 
discontinuity. Such two shocks appear in the winds from of the pulsar and stellar 
side [89]. The conditions at both shocks can differ significantly (magnetic field 
strength, acceleration efficiency) resulting in different acceleration efficiencies of 
particles and their energy loss rates. It has been shown that electrons accelerated 
on both shocks can reach different maximum energies saturating at ~10 GeV at 
the shock from the side of the massive star and up to several TeV at the shock 
from the pulsar side [89]. In summary, it is not surprising that precise predictions 
of the behaviour of the y-ray light curves from the massive binaries are very diffi- 
cult and they are not well described by the present models. The light curves and 
spectra may change in time due to not well known conditions within the binary 
systems discussed above. They depend on many parameters which at present are 
not constrained with enough precision. 

5.2.2. Hadronic models 

The production of y-rays in the massive binaries in the hadronic processes 
seems to be less likely due to a relatively long energy loss time scale of rela- 
tivistic hadrons in the hadronic collisions with the matter of the stellar wind (for 
the expected wind mass loss rate of the order of a few 10~ 7 M yr _1 ) with re- 
spect to their escape time scale with the plasma flowing from the binary sys- 
tem a long the shock struc ture. However, such models have also been considered 



(e.g. II253L 1254 Il92l I135IP . Definitive answer to the question of the importance 
of hadronic y-ray production in TeV binaries will come from observations with 
large scale neutrino telescopes by which the predicted neutrino rates might be 
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detectable (e.g. ni36ll288ll228lr >. 



5.3. Binaries with pulsars - prime targets for CTA 

Unexpected features of the y-ray emission from the TeV binaries have been 
already well documented by the present Cherenkov and satellite observations. 
The next generation of the Cherenkov arrays should be after all able to estab- 
lish whether the presently observed TeV binary systems are persistent sources of 
the modulated y-ray signal. If also TeV binaries other than LS I 61 303 show long 
time variability, what features determine the high level of the y-ray emission? Are 
these effects related in some way to the conditions created by the companion stars, 
e.g. similar to the long term activity observed in the Sun (e.g. the change of the 
stellar wind parameters - density, velocity, magnetic field - stellar cycles)? This 
problem could be answered by the multiwavelength observations in the high en- 
ergy domain together with the observations in the radio to optical range. A precise 
determination of the y-ray light curves and the measurement of the y-ray spectra 
at least in the Cherenkov energy range (sub- to multi-TeV) with an order of mag- 
nitude better sensitivity of the CTA will allow them for a more detailed testing of 
different cascade scenarios for the y-ray production mentioned in Sect. 5.2. 

CTA should also be able to measure shapes of the y-ray spectra at a few tens 
GeVs, i.e. in the region of the dip separating the high energy and very high energy 
components. Such detailed spectra at the sub-TeV energies should enable theo- 
retical studies of the origin of this feature (is it due to the presence of different 
components or other phenomena such as e.g., the absorption effects of y-rays in 
the UV-X-ray radiation ?). CTA will be able to make precise measurements of the 
phase dependent spectra from a much larger sample of the TeV binaries. Thus, it 
will allow us to investigate of the acceleration processes of particles at different 
conditions at the shock region. 

The TeV spectra, reported from the y-ray binaries, extend up to at least a few 
TeV without any evidence of the final cut-off. A detection of such cut-offs would 
allow us to constrain maximum energies of the accelerated particles. If these 
maximum energies extend up to a hundred TeV then the production of y-rays 
by hadrons is preferable since electrons in the environment of the binary system 
(strong magnetic and radiation fields) would have problems with overcoming the 
radiation energy losses. Therefore, the measurements of the high energy cut-offs 
in the y-ray spectra of these binaries with CTA would be of great importance for 
constraining the model parameters. 

The models of y-ray production in the pulsar winds predict the existence of 
a very hard y-ray component (spikes) being as a result of the comptonization of 
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the e xternal radiation field by leptons in the relativistic pulsar wind (e.g. [108l [77 



19511 . More recently, Cerutti et al. II 12511 compared the y-ray emission expected 



in the case of LSI 61 +303 and LS 5039 to the TeV observations of these two 
sources. It was concluded that the energies of e ± pairs in the pulsar winds cannot 
exceed 10-100 GeV. Future searches of such features in the spectra of TeV binaries 
(maybe only in some specific phases of the binary system?) with CTA (much 
better energy resolution and sensitivity in broader energy range) should enable 
unique determination of the basic parameters of the pulsar wind such as, e.g. its 
Lorentz factor, or the number of relativistic e pairs. 

Sensitive observations of the binary systems at sub-TeV energies with the next 
generation CTA might also make possible a detection of the pulsed component 
originating in the inner magnetospheres of the binary pulsars as observed in the 
case of the Crab pulsar. In the extremely compact binaries, containing a pulsar 
and the WR type star, this pulsed y-ray emission could be additionally modulated 
with the period of the binary system due to geometrical effects in the IC scattering 
process. Therefore, searches for new TeV binaries, containing energetic pulsars 
and WR type stars, such as that present in Cyg X-3, will be of great importance. A 
comparison of the emission features of different objects will enable a determina- 
tion of the crucial parameters determining the y-ray emission from these sources. 

Stellar and pulsar winds from y-ray binaries should produce large scale nebu- 
lae around the massive binary systems, filled with energetic particles accumulated 
during the lifetime of the pulsar. In addition to the observed y-ray emission mod- 
ulated with the period of the binary, such nebulae could also produce a steady low 
level y-ray emission llll5ll . Detection of this low level persistent TeV emission 
with CTA will be of great importance allowing us to constrain the history of the 
activity of this type of sources. 

It is difficult to reach clear conclusions on the number of y-ray emitting bi- 
naries within the Galaxy based only on the evolutionary studies due to unknown 
role of the energy output form the pulsar on the evolution of the companion star. 
Let us hope that the CTA will discover much more such objects during its scans 
of the galactic plane enabling population studies of the TeV binary systems in the 
future. 

In the case of extended binaries, of the LS2883/PSR1 259-63 type, a detection 
of the persistent and pulsed y-ray emission at phases farther from the periastron 
passage would allow us to investigate the direct role of the pulsar in the emission 
process when it is also close to the companion star. Therefore, CTA should try 
to measure sub-TeV emission from PSR1 259-63 when it is far away from the 
massive star. It will be interesting to find out whether the flaring GeV component, 
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discovered recently by Fermi-LAT about a month after the periastron passage, 
has persistent nature and similar features at different passages of the pulsar. A 
discovery of an accompanying sub-TeV component by CTA at these phases would 
make possible to constrain the particle spectra responsible for this component and 
put some limits on the emission region within the binary system (e.g. constrain 
the possible shock extension to regions very close to the companion star surface as 
predicted at some phases of the binary by Sierpowska-Bartosik & Bednarek [|267ll . 
This would provide new information about the nature of this surprising emission 
- whether it is due to an instantaneous state of the medium around the companion 
star or rather it is related to the properties of the pulsar and its wind, e.g. due to 
the geometrical effects in the collimation of the wind at specific directions? 



6. Microquasars 

The matter accreting onto a compact object such as neutron star or a black 
hole is expected to form an accretion disk. In the case of the accreting NS a part 
of the matter from the disk can be expelled forming a jet which can move with 
relativistic velocities. A similar jet is expected when the accretion occurs onto a 
rotating Black Hole (BH). Then, the rot ating BH i s expected to provide efficient 
mechanisms for energy extraction [|105[ l210l Il04ll . Jets are also expected to be 
accelerated and collimated hydromagnetically in the inner parts of the accretion 
disks [|208|] . An important consequence of these processes is the production of 
jets which moving relativistically in the direction perpendicular to the plane of 
the accretion disk. In fact, the large scale radio jets, some of them also showing 
featu res of the superlumi nal motion, were discovered in the galactic X-ray bina- 
ries 1122 U I222L 12851 Il85ll . Up to now, several such objects have been classified 
as microquasars in analogy to the quasar phenomenon observed around super- 
massive BHs in the nuclei of massive galaxies. Due to this analogy and also the 
evidence of the non-thermal radio to X-ray emission, indicating the presence of 
TeV electro ns (s ee e.g. observation of X-ray synchrotron radiation at parsec dis- 
tances II 13911 14011 . microquasars have been suspected to be sites of acceleration of 
particles able to produce y-rays in the GeV-TeV energy range. 



6.1. Observations of Gamma-rays from microquasars 

The possible detection of the y-ray emission from the well known binary sys- 
tem Cyg X-3, belonging to the mi croquasar class, were already claimed in data 



croqi 

from the past satellite telescopes 1 2041 122311 . but see [|179ll . Only recently, the 



GeV emission was reported in [|282[ 1 1 0f | . The observed y-ray emission is clearly 
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Figure 9: Fermi y-ray light curve (during the outburst) and the RXTE ASM X-ray light curve 
(entire lifetime of RXTE) observed from Cyg X-3 binary system ifioll . reproduced by permission 
of the AAAS). Note the anticorrelation between the GeV and X-ray emission. 



transient It appears close to the quenched radio state just be- 

fore the major radio flares observed in this system. The spectrum above 100 
MeV is well described by a single power law with the differential spectral index 
-2.7 ±0.05(stat) ±0.20(syst) and the peak flux ~ 2 x 10~ 6 photon cm~ 2 s _1 in 
the case of Fermi-LAT telescope | j 1 Of] and by the spectral index —2.02 ± 0.28 in 
the case of AGILE telescope [|l2lll . The emission shows clear modulation with the 
orbital period of the binary system, with the maximum close to the location of the 
compact object behind the companion star fl 1 Oil . This is also the case of the TeV 
binaries. The y-ray light curve is clearly anticorrelated with the X-ray light curve 
(see for details Fig. 9) Due to these similarities, it might be supposed that Cyg X-3 
should also emit TeV y-rays at phases when the compact object is in front of the 
companion star. Up to now, the TeV y-ray emission has not been discovered from 
Cyg X-3, in spite of the extensive observations with the MAGIC telescope 115 511 . 
However these observations start to constrain the y-ray spectra observed at low 
energies (see Fig. 10). 

In the case of another microquasar, within the massive binar y sy stem Cyg X- 
1, there is also evidence of a transient GeV-TeV emission [|52l |258Q. A flare of 
the TeV emission from Cyg X-l occurred close to one of the flares in the hard 
X-rays observed within the 3 day flaring episode (e.g. [121411 '). It lasted for less 



35 



■ Average Fermi/LAT active GbV period 4 soft state 

Fermi/LAT peak ^mission i i 

AGILE peak emission 

MAGIC soft state! ! 

MAGIC simultaneous GeV peak flux 



uL 



uL 



uL 



E [GeV] 



Figure 10: The y-ray spectrum from Cyg X-3 measured by Fermi (Abdo et al. 2009g), AGILE 
(from jl2lll and the upper limits from MAGIC B55|] , reproduced by permission of the AAS). Note 
that the upper limits from MAGIC start to constrain Cyg X-3 spectrum at TeV energies. 



than an hour when the compact object was behind the massive companion. The 
significance of this flare is rather low (below 5 a). Also a separate flare abov e 10 



MeV, on a time scale of 1-2 days, was reported recently by the AGILE Team [258]. 
The GeV flares seem to be rather exceptional during the 300 day observation 
period with the AGILE telescope Ill47ll . Unfortunately, the GeV flare has not been 
detected by Fermi-LAT telescope during this same period. The y-ray emission 
from microquasars in the low mass binary systems have not been discovered up 



to now in s pite of ex tensive observations, see e.g. GRS 1915+105 [29. 125911 or 



SS433 [136|,ll78, 259]. 



6.2. Theoretical modelling and expectations 

It is generally expected that the high energy y-ray emission in the previously 
discussed two microquasars is related to the presence of a massive star close to 
a mildly relativistic jet. A few different sites within the massive binary systems 
have been discussed as responsible for the origin of y-ray s. Close to the base of 
the jet, where the magnetic field strength is expected to be large, nuclei are pref- 
erentially accelerated to energies enabling the GeV- TeV y-ray production. How- 
ever, electron acceleration is expected to be saturated there at relatively low values 
(GeV range) due to the synchrotron energy losses, thus preventing the TeV y-ray 
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production. On the other hand, the acceleration of nuclei should not be strongly 
influenced by their energy losses allowing them to reach even PeV energies. Nu- 
clei with such energies may interact efficiently with the X-ray pho tons of the inner 
parts of the accretion disk producing y-rays and neutrinos (e.g. Q207I]). At inter- 
mediate distances from the accretion disk, its radiation field is still strong enough 
for an efficient disintegration of nuclei but it is too weak for the pion production in 
collisions with photons [82]. In this case, neutrons released from the nuclei, can 
outshine the accretion disk producing y-rays and neutrinos in collisions with the 
matter of the disk (see e.g. I1227ll ). At larger distances from the disk, the radiation 
field of the companion star starts to dominate over the radiation field of the ac- 
cretion disk. Then, neutrons from the disintegration of nuclei outshine mainly the 
massive companion producing neutrinos and y-rays in collisions with the stellar 
atmosphere [82]. At these distances the synchrotron energy losses of accelerated 
electrons become low enough, allowing them to be accelerated to TeV energies. 
Therefore, electrons can efficiently prod uce y-rays by the inverse Compton scat- 
tering of the stellar radiation (e.g. [|l67[ 252 1, and recently | 311 ]. The outer re- 
gions of the jet (which are still within the binary system) seem to be likely sites 
for the production of GeV-TeV emission which fulfills the observational charac- 
teristics obtained from the observation of the Cyg X-3 binary system, i.e. the 
modulation of the GeV signal with the period of the binary system. However, 
y-ray production in the interaction of particles with the radiation and matter of the 
accretion disk can also be an important mechanism, provided that the observed 
low level of modulation of the y-ray signal is due to other effects, e.g. the preces- 
sion of the disk and/or the jet with the period of the binary system. 

The high energy y-rays produced in the compact binary systems, containing 
massive stars (such as Cyg X-l and Cyg X-3), sho uld a lso suffer from the strong 
absorption in the radiation field of the star (e.g. [|224ll ). As a result, the spec- 
tra of the escaping y-rays are formed in the IC e ± pair cascade, developing in 
the anisotropic radiation field of the massive companion (as in the case of TeV 
binaries, see Sect. 5.2.1). The detailed calculations of the optical depths of the 
massive star environments in the binary systems Cyg X-l and Cyg X-3, show that 
the escape of TeV y-rays depends strongly on the phase of the binary system, see 
e.g. [92] for Cyg X-l and [|88|l for Cyg X-3. It was shown that the strong TeV 
y-ray emission from Cyg X-l at the phase reported by the MAGIC observations 
i.e. at the location when the compact object is behind the massive star 1I52I1 . is 
rather unlikely (19211 . T he in jection spectrum of electrons in this case should be 
very hard, as shown in [130911 . Due to the very strong y-ray absorption close to the 
hot WR type stars, it is also rather unlikely to expect the present Cherenkov Tele- 
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scopes to detect any emission from Cyg X-3 [88]. In principle, the y-ray emission 
observed from Cyg X-3 might originate also as a result of the ICS of the disk 
radiation by relativistic electrons in the inner jet [U27I1 . Detailed calculations of 
the y-ray spectra along such general scenario, which include the development of 
the anisotropic IC e ± pair cascade in the disk radiation, has been recently made 
by Sitarek & Bednarek 11270 1. It is concluded that the observed y-ray spectrum 
from Cyg X-3 can be explained provided that the emission region in the jet is at 
the distance up to 100 — 300 inner radii from the accretion disk. In this case, the 
observer should be located at a relatively small angle to the jet axis. 

Relativistic particles may appear in the binary system due to acceleration on 
the shock waves being a result o f the interactio n of the jet with a dense, clumpy 
wind from the massive star (e.g. [|240ll255[|l55ll . In such scenario y-ray s might be 
produced in hadronic collisions, as considered by e.g. 11691 123311 . The y-ray fluxes 
strongly depend on the properties of the stellar wind (density, velocity) which are 
at present not well known. Also the interaction of the jet in its termination region 
(on a parsec scale) may result in the production of a steady , high e nergy y-ray 
emission with an intensity difficult to determine reliably, e.g. [|73Llll2|l . However, 
a discovery of such a persistent GeV-TeV y-ray emission would provide important 
constraints on the content and energetics of the microquasar jets. 

Originally, the TeV y-ray binaries (see Chapter. 5) have been considered to 
belong to the microquasar class and their early modelling postulated the injection 
of relativistic particles along the jet (e.g. 1I1491 l85ll . However, now the hypothesis 
of the presence of energetic pulsars in these systems is more widely accepted. 



6.3. Open questions to be answered by CTA 

Interesting questions concerning the acceleration of particles, their propaga- 
tion and radiation process in microquasars may be put in new light with an order 
of magnitude more sensitive observations above a few tens of GeV. Long term 
monitoring at sub-TeV energies should confirm the observed relation of the y-ray 
emission only to the major radio flares of Cyg X-3. However, the question also 
appears whether lower level y-ray emission is produced also in the quiescent ra- 
dio stages? In the case of the high mass microquasars, such as Cyg X-3 and Cyg 
X-l, the localization of the emission region within the binary system is of cru- 
cial importance. An observation of a break in the y-ray spectrum at a few tens of 
GeV would be consistent with the absorption by the stellar radiation in the sys- 
tem. However, the predicted fluxes on the level of ~ 1% CU at 20-30 GeV and 
0.1% CU above 100 GeV §£M would be probably difficult to detect even with 
CTA particularly that the emission has a transient nature with the duty time of the 
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order of a few days. A cut-off in the spectrum at even lower energies might sug- 
gest y-ray absorption by the radiation of the hot inner accretion disk, indicating 
the localization of the emission region in the inner jet. It will be interesting to 
search with CTA for an anticorrelation of GeV and (possible) sub-TeV fluxes as 
observed in the case of TeV binaries. Their discovery would show an importance 
of the y-ray absorption and a formation of "escaping" y-ray spectra in the cascad- 
ing processes. A confirmation of the detection of the TeV emission from Cyg X-l 
would be of great importance since stars in Cyg X-l binary system have much 
larger separation than those in Cyg X-3. Thus, the escape conditions for the TeV 
y-rays produced in the binary system are much less restrictive. 

CTA should be able to conclude whether microquasars in the low mass X-ray 
binaries (LMXBs) can also produce y-rays. The acceleration of particles to TeV 
energies in LMXBs is expected to occur also efficiently as conclud ed from th e 
observations of synchrotron radiation in the X-ray energy range (e.g. ||139lI140I1 ). 
This process is expected to be determined rather by the jet properties than by 
the companion star. The upper limits (or low level detections) of the GeV-TeV 
photons with CTA would allow us to answer the question on the importance of 
a massive star in the microquasar for high energy processes. For example, in the 
LMXBs y-rays might be produced in another mechanism than that expected in the 
HMXBs, e.g. the Synchrotron Self-Compton process might dominate. 



7. Massive stars in the close binary systems 

Strong winds of massive stars (O and WR type) may collide producing a dou- 
ble shock structure separated by a contact discontinuity. These winds are charac- 
terised by large energies, large velocities and strong magnetic fields (due to the 
strong surface magnetic fields of the companion stars). It is expected that these 
shocks accelerate particles (electrons and hadrons) to energies allowing them to 
produce GeV- TeV y-rays. These systems seem to be unique laboratories for the 
investigation of the mechanisms of particle acceleration (shock acceleration or re- 
connection of magnetic fields) due to the well defined target for the relativistic 
particles within the binary system (density of matter, soft radiation field). 

7.1. High energy observations of massive binaries 

Up to now only one massive binary system, Eta Carinae, has been found in th e 
direction of the Fermi-LAT source at energies below ~100 GeV ll282l I20L 161 1. 
The high energy emission of this source is characterised by an intriguing two- 
component spectrum which can be approximated by two power laws, the first one 
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Figure 11: The y-ray spectrum observed from the direction of Eta Carinae binary system. Two 
spectral components are clearly observed, the soft o ne, s howing exponential cut-off at a few GeV 
and the hard one, extending up to ^100 GeV (from 016211 . reproduced by permission of the A&A). 



with an exponential cut-off at a few GeV and the second extending up to at least 
~100 GeV (see Fig. ITTb . Recently, it was shown that the high energy part of 
this emission must be associated with the Eta Carinae lll62ll . Whether this second 
component extends through the TeV energy range or cuts off earlier is at present 
unknown. 

It is also possible that a part of the TeV emission from the direction of the open 
cluster Westerlund 2 might originate in the massive binary system, WR 20a [45]. 
However, it is not yet clear whether any emission comes from the Westerlund 2 
or from the nearby PWNa D25J. Such PWNa is probably associated with the y-ray 
pulsar recently detected by the Fermi-LAT [|263ll . The upper limits on the TeV 
emission from other massive binary systems, containing WR type stars such as 
WR 146 and WR 147, have been reported recently by the MAGIC Collabora- 
tion [61]. They are on the level below the theoretical predicti ons a llowing one to 
constrain the parameter space of some models for their origin Q249Q. 



7.2. Models for y-ray emission 

The high energy radiation, expected in te rms of the wi nd co l lision mo del, was 



considered in more detail in several papers 

[HlioilElHHtMS]. In 

general, the acceleration process of electrons and hadrons can be saturated by 
the radiative energy losses or by particle escape from the binary system. In the 
case of electrons, the synchrotron process and the Inverse Compton scattering on 
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the stellar radiation play the dominant role. Electron acceleration can be satu- 
rated by the synchrotron losses at TeV energies. However, electron acceleration 
can be also saturated by IC energy losses at much lower maximum energies, of 
the order of several GeV (for detailed formulae see e.g. V96n ). The above men- 
tioned values are typical for binary systems containing two massive stars. Which 
of these saturation scenarios is realised in particular binary system, depends on 
the magnetic field and the radiation field of the massive star. On the other hand, 
acceleration of hadrons can be saturated either by energy losses in hadronic col- 
lisions with the wind matter or due to the escape from the binary system along 
the shock. The maximum energies of hadrons determined by hadronic interac- 
tions were estimated as a few hundred TeV (see, e.g. D96|]). On the other hand, 
the maximum energies of hadrons, if their escape from the binary system is im- 
portant, are expected to be of the order of a few PeV. It is clear that depending on 
the conditions at the shock, the acceleration of electrons and/or hadrons is satu- 
rated at very different maximum energies. This may lead to a variety of radiation 
scenarios depending on the phase of the binary system. They can be significantly 
different depending on the conditions in the binary systems. In order to put some 
light on the possible high energy processes occurring in the binary systems and 
on the acceleration mechanism of particles, only long term sensitive observations 
of these systems at high energies, simultaneous with those at other parts of the 
electromagnetic spectrum, are required. 

A few binary systems have been considered in more detail as possible y-ray 
emitters (e.g. WR 140, WR 146, WR 147, WR 20a, Eta Carinae). In general, 
the models postulating the acceleration of electrons do not predict y-ray emission 
at energies clearly above ~100 GeV due to a sa turation of the electron acceler- 
ation by radiative energy losses (e.g. [ 249 , 243fl). On the other hand, hadronic 
mod els pr edic t y-ray e mission which in principle can extend up to PeV ener- 
gies [| 102L l8ll ll62L l96|l . However, this emission, if produced in the binary sys- 
tem, should be st rongly influenced by the absorption in the thermal radiation of 
the stars 11821, l249i 19611 . Moreover, y-ray s with TeV energies can be produced also 
by hadrons escaping from the binary system into the relatively dense nebulae sur- 
rounding it (e.g. 1123 ill ). This emission should be steady in contrast to the y-ray 
emission coming from the binary system which is expected to be modulated with 
the binary period due to the eccentric orbits of the stars and geometry dependent 
absorption of y-ray s with energies above ~100 GeV. Therefore, various y-ray 
emission features are expected from massive binary systems. 
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7.3. Future studies with CTA 

A variety of the emission features predicted by different scenarios can only be 
tested by more sensitive instruments such as CTA. CTA will be able to provide 
spectral information in the whole range of the binary phases i.e. at different con- 
ditions in the acceleration region. No detection of any y-ray signal above ~100 
GeV would support their leptonic origin in the binary systems. However, a detec- 
tion of y-rays above ~100 GeV, modulated with the binary period, would support 
the acceleration of hadrons within the binary. 

It would be interesting to determine whether the acceleration of electrons 
and/or hadrons occurs only at specific phases, depending on the conditions at the 
shock. Such investigations should throw new light on the conditions for acceler- 
ation of electrons or hadrons at shock waves. On the other hand, CTA could also 
discover a weaker steady y-ray component from the surroundings of the binary 
systems containing two massive stars, thus enabling studies of the particle diffu- 
sion. This will put additional constraints on the acceleration process of particles 
at the shocks with different conditions. 

8. Accreting White Dwarfs 

Accreting magnetized White Dwarfs can provide a variety of conditions for 
acceleration of particles either in the inner, turbulent part of the accretion disk 
(the accretor or propeller phase) or in the thermonuclear explosions occurring on 
the surface of the WDs. Thermonuclear explosions may occur from time to time 
in the matter accumulated on the WD surface. The expelled matter creates a shock 
able to accelerate particles. 

Fermi-LAT has recently discovered y-ray emission from the symbiotic binary, 
V407 Cyg, containing a WD orbiting a red giant. The y-ray flare is correlated with 
the optical Nova tl9fl. The observed spectrum is hard below ~1 GeV showing 
an exponential cut-off at larger energies (see Fig. [12]). An extrapolation of the 
spectrum (from the two last points with a simple power law) is above the 5 hr 
sensitivity of the planned CTA at energies above ~100 GeV [30]. 

The mechanism of this y-ray production (leptonic or hadronic) is at present 
unknown. The observed y-rays are believed to be produced either by decay of 
pions produced by hadrons interacting with a dense wind of the red giant or by 
electrons which IC scatter the soft photons from the red giant. If it is leptonic, 
then additional high energy component in the spectrum could be present due to the 
accelerated hadrons. We conclude that future Nova explosions, although expected 
to be rather rare (the rate between 0.5 — 5 yr 1 in the Galaxy 121111 . have a chance 
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Figure 12: y-ray spectrum from the direction of the symbiotic binary V407 Cyg containing a White 
Dwarf on an orbit around the red giant (from 11191 . reproduced by permission of the AAAS). Note 
that the spectrum seems to cut-off exponentially above ~1 GeV. 



to be detected by CTA. Such detections will put important constraints on the shock 
acceleration model at the conditions met in the symbiotic binaries and on the y-ray 
radiation mechanisms. 

Recently, y-ray emission from binary systems containing accreting WDs was 
studied in j219L I95Q. In these models, it is proposed that the turbulent transition 
region between the rotating magneto sphere of the WD and the accreting matter 
can provide good conditions for the acceleration of particles (the propeller stage). 
These particles might be able to produce even TeV y-rays in some systems, ex- 
pected to be detectable by Cherenkov telescopes. However, accelerated electrons 
are probably not able to produce y-rays with energies above ~10 GeV due to the 
saturation of their acceleration by the synchrotron radiation. 



9. Other objects expected to emit TeV y-rays 

The high energy processes in other Galactic objects have also been analysed 
in terms of the leptonic and hadronic models. These objects, e.g. isolated massive 
stars, White Dwarfs or Young Stellar Objects and accreting compact neutron stars 
(LMXBs) and White Dwarfs, have been predicted to emit the y-ray fluxes on the 
level detectable by the future Cherenkov telescopes. 
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9.1. Winds of the classical massive stars 

Strong winds, produced by the early type stars (Be, O, WR), are expected to 
be non- stationary, containing smaller scale shocks. It has been argued that these 
shocks can accelerate particles ei ther close to the stars rl3Q4L l244h or at the terminal 



shocks of the stellar winds (e.g. BOOl 112811 ). Electrons accelerated in such shocks 



have energies allowing them to produce of GeV y-rays in the bremsstrahlung pro- 
cess Il244n . On the other hand, had rons accelerated to ~100 GeV could produce 



y-rays via decay of pions (e.g. 1130511 . The predicted fluxes seemed to be within the 
sensitivity of the EGRET telescope but no source of this type had been reported 
at that time. Massive stars usually stay close to their birth places which are within 
massive stellar associations. They are also expected to be sources of y-rays due 
the presence of large scale shocks produced by a accumulation of stellar winds. 

In fact, the TeV emission was detected from the directions of a few open stellar 
clusters: Cyg OB2 [Ell], Westerlund 2 0,0 and Westerlund 1 ll2~3~oll . Its origin 
is still debated. It is not excluded that a part of this emission comes from isolated 
massive stars. However, it will be rather difficult to identify the TeV emission 
from specific massive stars within the open clusters even with the CTA angular 
resolution and sensitivity due to a large number of potential sources (several mas- 
sive WR type stars, pulsars, supernova remnants, etc ...). However, O and WR 
type stars, which are sometimes well outside clusters, can become interesting tar- 
gets for the observations with CTA. Their possible detection will enable studies of 
the high energy processes occurring in the winds not far from stellar surfaces, con- 
straining the structure and strength of the wind magnetic fields and high energy 
flaring activity in the magnetospheres of massive stars. 

9.2. Young stellar objects 

An other promising candidates are the massive young stellar objects (YSO) 
which create asymmetric outflows (jets) propagating up to the distances of a frac- 



tion of pc R217I1 . These jets may collide with dense molecular clouds forming bow 
shocks. If a YSO is located within a molecular cloud (e.g. IRAS 16547-4247 
or HH80), then hadrons, interacting wit h the distributed matter, could produce 



TeV y-ray fluxes observable by CTA [|68[ 111411. T he non-thermal radio emission, 



detected from the field of YSO BD+433654 || 10311 . indicate the importance of non- 
thermal processes in these type of sources. A modelling of such emission predicts 
the production of TeV y-rays above the sensitivity limit of CTA. 

A transient acceleration of hadrons to TeV energies is also expected to occur in 
the magnetospheres of the accretion disks of low mass young stellar objects such 
as T Tauri stars fl 14811 . These particles, colliding with matter, may produce y-rays 
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which, after absorption in the stellar and disk radiation, can escape only with en- 
ergies below ~100 GeV. Although, the expected y-ray luminosities are relatively 
low, T Tauri stars can be promissing targets for the next generation Cherenkov 
telescopes. The transient nature of their y-ray emission may however lower the 
detection probability. Note that there seems to exist a correlation between the 
y-ray objects in the Fermi-LAT catalogue and those in the catalogue of known 



YSO 122511 . 

Efforts to study massive young stellar objects and the T Tauri stars should 
be certainly undertaken with the sensitivity of CTA. These objects represent dif- 
ferent conditions for possible acceleration of particles (e.g. low velocity jets but 
relatively dense medium). Their possible detection will provide interesting con- 
straints on the acceleration processes of particles and formation of jets in con- 
ditions significantly different from those expected in the sources with relativistic 
jets. 

9.3. Isolated White Dwarfs 

The processes for the high energy y-ray production in the inner magneto- 
spheres of rotating neutron stars (discussed in Sect. 2) may in principle operate 
also in the vicinity of the fast rotating White Dwarfs (WDs). Strongly magne- 
tized, fast rotating WDs are expected to be a final stage of the evolution of the 
Ap and Bp type stars or as a result of mergers of WDs in the binary systems. In 
fact, strong magnetic fields have been discovered in the case o f a fe w percent of 



isolated WDs and ~25% of WDs in the Cataclysmic Variables [306]. 

The possibility that magnetized WDs can lose large amounts of the rotational 
energy on the generation of electromagnetic waves and acceleration o f particles 



was considered i n the cont ext of two binary systems: IE 2259+586 112341 . 12921 



and AE Aquarii 111421, 118911 . It is expected that magnetized WDs are able to in- 
ject relativistic electrons with energies allowing them to produce TeV y-rays in 
the IC process (see e.g. IU91I1 '). In fact, a transient TeV emission from the binary 



system AE Aqr, containing fast rotating WD, was reported in th e past II 11 61 . 1 1291 



These reports were not confirmed by more recent observations [205. 126511 . How- 
ever, due to its likely transient nature, the lack of confirmation does not exclude 
possible emission of y-rays by such type of objects. Systematic observations of 
magnetized WDs with arrays of the next generation Cherenkov telescope might 
result in the detection of magneto spheric y-ray emission from WDs. These y- 
rays could be produced as a result of ICS of thermal radiation from the hot stellar 
surfaces of the WDs. Note that due to the significantly weaker magnetic field in 
the inner magnetospheres of the rotating WDs, such pulsed y-ray emission should 
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extend to larger energies (GeV-TeV range) than those observed in y-ray pulsars. 
Moreover, steady TeV y-ray emission, produced by leptons escaping from inner 
magnetospheres, could also appear as a result of the comptonization of the CMB 
or the infrared to optical photons produced in the Globular Clusters containing 



hundreds to thousands of magnetized WDs [90]. 

The sensitivity of CTA gives a chance to detect y-rays from the nearby fast 
rotating magnetized WDs at distances up to ~ 300 pc (assuming typical period of 
100 s, the surface magnetic field of 10 8 G and the efficiency for y-ray emission of 
10%) H- 

9.4. Accreting Neutron Stars (LMXBs) 

Old NSs in binary systems have usually lost most of their rotational energy. 
Due to the long rotational periods, they are not able to create winds strong enough 
to push out the matter outside their light cylinder radii. They accrete matter from 
the companion star in the two basic regimes, the so-called propeller and accre- 
tor stages. In the propeller stage, the compact object accretes matter only up to 
a certain distance from its surface, where the magnetic pressure of the rotating 
magnetosphere balances the pressure of the accreting matter. Such scenario was 
discussed in the past in the cont ext of possible high energy processes in the so- 
called hidden pulsar model [ 279 . 281 1 and in the cauldron model ll99l . l290ll . More 



recently y-ray emission from binary systems containing accreting NSs was con- 



sidered in Q86J]. It was concluded that the turbulent transition region between the 
rotating magnetosphere of the compact object (propeller stage) and the accreting 
matter could provide good conditions for the acceleration of particles. In some 
systems, the particles might be able to produce even TeV y-rays detectable by the 
Cherenkov telescopes. 

In the case of the accreting objects in compact Low Mass X-ray Binaries 
(LMXBs), huge gravitational energy of accreting matter is converted mainly into 
X-ray emission. It has been suspected that some of this power might also go to 
y-rays. Therefore, LMXBs have been considered as potential targets for y-ray 
telescopes. X-rays produced close to the accreting object or in the accretion disk 
can heat the nearby companion star to a temperature significantly above that ob- 
tained from nuclear burning. Therefore, it is likely that the irradiated stars produce 
soft radiation fields for relativistic electrons comparable to those observed in the 
massive TeV y-ray binaries (see [94J. Up to now, none of the LMXBs have been 



detected by the modern Cherenkov telescopes (see e.g. [56]). However, the first 



object belonging to the class of compact Low Mass X-ray Binaries might be al 



ready detected in the GeV energies by the Fermi-LAT telescope [274]. These 
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objects remain interesting targets for the future instruments with good sensitivity 
above a few tens GeV such as CTA. 



9.5. Binary systems containing MSPs 

A significant part of the millisecond pulsars (MSPs) are also found in the com- 
pact binary systems. As discussed above for classical pulsars, also MSPs can pro- 
duce winds which interact with winds of compan ion stars p roducing shocks able 



to accelerate electrons to GeV-TeV energies (e.g. [I197U242L I7 ill. The y-ray fluxes 



from such binary systems are expected to be low due to much weaker radiation 
fields of the companion stars in these systems. The soft radiation field from the 
surface of the companion star can be larger than that produced in nuclear burning 
in the case of a strong illumination of the stellar surface by X-rays from the pulsar. 
Then, the companion star may be heated to significantly larger temperatures than 
those expected from the pure nuclear burning mechanism (e.g. [94]). Relativis- 
tic electrons from MSPs, interacting with this enhanced stellar radiation, would 
produce y-rays as a result of IC scattering. Detailed y-ray light curves from such 
binaries are expected to differ significantly from those observed in the TeV y-ray 
binaries due to the anisotropic emission of thermal radiation by the companion 
star. 

Some of the binary systems containing MSPs can be completely surrounded 
by the stellar winds creating the so called "hidden" MSPs ll279Tl . In principle, such 
MSP binaries can also contribute to the y-ray emission observed from individual 
objects or from Globular Clusters. 



10. Prospects for detection/identification of the Galactic y-ray sources with 
CTA 

About a few hundred GeV y-ray sources in the Galactic plane [? ] and several 
TeV y-ray sources (see TeVCat Catalogue and also [48 J) have no clear associa- 
tions with known objects in the Galaxy. Many of these sources observed in GeV 
energies are probably pulsars and in the TeV energies the PWNe surrounding 
them. In fact, out of the 25 PWNe in the TeVCat Catalogue, 16 have associations 
with the 2nd Fermi Catalogue [? ]. Many of the TeV unidentified sources show 
some angular extension what strengthen the hypothesis of their relation to PWNe 
or SNRs. On the other hand, out of the unidentified TeV sources, 17 are associated 
with the 2nd Fermi Catalogue sources. The lack of clear objects, seen in radio or 
X-rays, in the directions of these TeV sources puts strong constraints on the y-ray 
radiation mechanism (leptonic, hadronic ?), and also on the transport process of 
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particles in the sources and on their evolution. Therefore, detailed investigation of 
this broad class of unidentified GeV-TeV y-ray objects, with an order of magni- 
tude more sensitive instrument and with the better angular resolution (CTA), will 
provide important constraints on the sites of particle acceleration in the Galaxy. 

Significant amount of the unassociated GeV sources tend to group in the direc- 
tions to regions in the Galaxy with bright diffuse y-ray emission. These sources 
often show curved spectra (about 28% [? ]) suggesting that their TeV emission 
may be on a much lower level. One of the main tasks of CTA will be to pro- 
vide information on the shape of the y-ray spectra at sub-TeV energies of these 
sources in order to conclude on the spectral behaviour: Do they really cut-off or 
do they show a second high energy spectral component? For this purpose, the fu- 
ture instrument should have a low energy threshold in order to obtain good quality 
spectra at energies clearly below ~100 GeV. The CTA superior sensitivity and 
good angular resolution should enable finding associations of these y-ray sources 
with the objects known in other parts of electromagnetic spectrum or with already 
known classes of the y-ray sources (e.g. SNRs, PWNe, binaries, etc ...). 

The large field of view of the CTA would make galactic plane scans possible 
with an order of magnitude better sensitivity than it has been done up to now by 
the HESS collaboration J37 L 43] and also in the case of the Cygnus region by the 
VERITAS Collaboration Il63ll . Such scans should be done regularly in the same 
parts of the Galactic plane in order to have a chance to discover variable and pe- 
riodic y-ray sources. Up to now, a number of such variable known sources is 
relatively small probably due to the limited sensitivity of the present Cherenkov 
telescopes. Frequent scanning of the same regions of the Galactic disk with CTA 
should result in discovery of new classes of variable TeV y-ray sources not ob- 
served up to now. 

CTA should also have the power to investigate fine structures of the much 
extended sources of the type recently reported by the MILAGRO detector (e.g. 
MGRO J2019+37, MGRO J2031+41, C3, or C4 ELI!]). An investigation of the 
specific parts of these sources (spectral information, possible variability ?) would 
allow us to conclude whether these TeV sources are due to either the composition 
of a few nearby different types of sources (e.g. in the Cygnus region) or they rep- 
resent separate parts of a nearby source (e.g. source C3 coincident with Geminga 
pulsar?). With the present ARGO-YBJ and the planned construction of the next 
generation of the water tank Cherenkov detectors (HAWC and LHAASO), the 
number of such extended sources will certainly increase being a very interesting 
class of objects to be studied with CTA. 
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11. Conclusion 

The important role of the high energy processes in a few types of compact 
Galactic sources have been recently undoubtedly established with the discovery 
of a few hundred GeV and a few tens of TeV y-ray sources. Many surprising 
emission features have been already discovered. Their detailed investigation will 
certainly put a new light on the origin of the high energy particles in the Galaxy. 
The nature of the high energy processes resulting in the production of the high 
energy y-rays remains mostly unknown and is at present debated. The appearance 
of the y-ray emission in the case of some sources have not been expected by theo- 
retical modelling. In order to understand these high energy phenomena, extensive 
deep studies with the next generation of the TeV telescopes are required. In fact, 
the scientific community realized these needs initiating new instrumental projects 
in the high energy y-ray astronomy such as CTA, HAWC, or LHAASO. They 
will be able to study the already discovered classes of objects with an order of 
magnitude better sensitivity. The main tasks of the future TeV observations will 
be: 

• Enlarge the number of sources in specific classes of y-ray objects in order 
to enable comparative studies. 

• Discover new types of sources. The importance of high energy processes, 
leading to the y-ray production, is still expected in specific objects based on 
their theoretical modelling. 

• Perform multiwavelength observations of the known types of sources (mainly 
in radio, X-ray, GeV and TeV y-rays) in order to get complementary infor- 
mation on their high energy features. 

• Perform morphological studies of the extended sources through the sub-TeV 
(~ 30 GeV) to multi-TeV (>30 TeV) energies in order to determine the sites 
of particle acceleration, their propagation conditions within the sources and 
their surroundings, and answer the question whether the specific parts of the 
extended sources have common origin (single source of relativistic particles 
and/or radiation mechanism) or appear by chance as a result of overlapping 
of different (types of) sources. For example, in the case of the open clus- 
ters: pulsars, supernova remnants, binary systems and microquasars can 
contribute simultaneously to the observed y-ray flux. 
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• Try to determine the highest energy end in the TeV y-ray spectra of the 
sources in order to get constraints on the maximum energies of accelerated 
particles and on the acceleration process. 

• Establish the shortest possible variability time scale of the TeV y-ray emis- 
sion in known (and new) types of sources in order to obtain constraints on 
the acceleration and cooling time scales of relativistic particles and perhaps 
answer the question on their production mechanism (hadronic or leptonic 
?)• 

• Study the spectral properties of the periodic y-ray sources in order to inves- 
tigate the acceleration of particles and their radiation processes at different 
conditions at the source (magnetic field strength, density of matter and ra- 
diation fields). 

• Study the long term variability of periodic sources (pulsars, binary systems) 
and try to link these features to the conditions at the source derived from 
other observations (e.g. variability in the radio light curve of pulsars or 
their giant flares, or stellar cycles (similar to Solar cycles) resulting in the 
change of stellar wind properties?). 

• Establish a relation between the y-ray sources observed in the GeV energy 
range (by Fermi-LAT and AGILE) to their properties in the TeV energy 
range (investigate spectral and morphological features). This should make 
possible to put constraints on the radiation processes engaged in the spe- 
cific sources (importance of the absorption at the source and cascading pro- 
cesses?) and determine the main production process (hadronic versus lep- 
tonic?). 

• Study a relation of the variable y-ray components in the sources (produced 
by particles accelerated recently) to the possible lower level steady compo- 
nents (produced in the source or around it by particles accelerated in the 
past) in order to constrain evolutionary properties of these sources. 

The outcome of the investigations summarized above will provide conclusive 
information on the physics of sources discussed here (such as the particle acceler- 
ation processes, properties of relativistic magnetized plasma, structure of rotating 
magnetospheres, formation of relativistic jets, propagation of particles, etc ...). It 
is hoped that, with the next generation of instruments, the basic question on the 
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origin of energetic particles in the Galaxy will be finally answered. At least the 
author of this paper is convinced that this will happen. 
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